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BENDING OF A RECTANGULAR PLATE WITH ONE FREE EDGE 


W. Nachbarl 
(Proc. Paper 1196) 


ABSTRACT 


A rectangular plate, simply supported along three edges and free on the 
fourth, is subjected on one face to a lateral pressure which is uniform paral- 
lel to the free edge but can vary linearly in a direction perpendicular to the 
free edge. Infinite series for deflections, moments, shear and reaction force 
distribution were derived and programmed for computation on a high-speed 
digital computer. Numerical values are presented in graphs and tables for 
the two special cases of completely uniform and triangular loading; the more 
general loading is given by superposition. The discussion of results includes 
(a) comparison of results for plates with different values for Poisson’s ratio, 
(b) computational difficulties and methods used to overcome them, and (c) 
approximate formulas in closed form. 


INTRODUCTION 


This note considers bending, according to the classical thin plate theory of 
Poisson-Kirchoff, as presented in Timoshenko’s standard text on plates, | 1! * 
of a homogeneous isotropic, rectangular plate which is free on one edge and 
simply supported on the remaining three edges. The contour of the lateral 
pressure loading on the plate is trapezoidal parallel to the simply supported 
edges and uniform parallel to the other pair (See Figure 1). Formulas in the 
form of infinite series are derived for deflections, bending and twisting mo- 
ment distribution, shear distribution, support reaction distribution and cor- 
ner reactions. These series are presented in a form which was found to be 
particularly convenient for numerical computation using a high-speed digital 
computer. 

The bending of an isotropic rectangular plate which has three simply sup- 

' ported edges and one free edge is treated briefly in Timoshenko (1) pages 

218-219, and a table of values is given there for deflection and maximum 


Note: Discussion open until September 1.1957. Paper 1196 is part of the copyrighted 

Journal of the Engineering Mechanics Division of the American Society of Civil 

Engineers, Vol. 83, No. EM 2, April, 1957. 

1. Research Branch, Missile Systems Div., Lockheed Aircraft Corp., Palo 
Alto, Calif. 

* Numbers in square brackets refer to corresponding items in the list of 

references below. 
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Figure 1. Geometry Loading and Edge Conditions for Plate. 
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bending moment under a uniform load. More recently, Goriupp [2] derived 
the Green’s Function for this plate and integrated this to obiain the solutions 
in series form for the plate under uniformly and under triangularly distri- 
buted lateral pressure. A numerical table of moment, shear and reaction 
values was presented in this paper only for a ratio of plate length to breadth, 
b/a, equal to 0.75. The main emphasis of this paper, however, was on the 
development of simple finite approximating formulas rather than on presen- 
tation of numerical results. 

The infinite series found by Goriupp were not used for present purposes 
for two principal reasons. First, they were derived by setting Poisson’s 
ratio “ equal to zero in the free edge boundary conditions; hence, the series 
were exact only for the case w= 0. The finite formulas which were developed 
involved the same assumption. Secondly, the series coefficient were defined 
using the hyperbolic functions. While such definition makes for compact 
notation, it leads to difficulty when computing such slowly converging series 
as for the edge reactions. The coefficients of high order in this case are 
discovered to be computed as the very small differences between very large 
numbers, with a consequent loss of numerical accuracy. This difficulty, and 
the step taken to avoid this for the present work, are further discussed in the 
following section. 

The analysis derives a single series solution for the general trapezoidal 
loading, the contour of which is governed by a parameter k (Equation 1). Us- 
ing the solutions for the two special cases of uniform load, k = 0, and trian- 
gular,load, k = 1, the solution for any general trapezoidal loading may be 
found by linear superposition. For the value & = 0.3, the derived series were 
computed for both loading cases and over a wide range of values for b/a; the 
results are presented graphically on a non-dimensional scale in Figures 8 
through 22 below. In addition, for the single value b/a = 0.75, the series were 
computed for several values of w including 0 and 1/2. The results are pre- 
sented graphically in Figures 3 through 7, and in Table 5, and are used to ob- 
serve the effect of Poisson’s ratio upon deflection and stress distribution. 


Formal Solution 


Let (x, y) be a cartesian coordinate system in the middle plane of the plate, ies 
and let w(x,y) be the lateral deflection* of the middle surface (Fig. 1). The 
lateral loading q(x,y) on the rectangular plate of side lengths a, b is assumed 
to be in the form a, hs 
a(x, (1) 
The differential equation for w(x,y) is then 
ow ° > (2) 
a ax* ay* oy 
For simple support along edges x = + a/2 we seek to represent w(x,y) by 4 
* The coordinate system is assumed to be right-handed with w and q positive ed ° 
in the + 2 direction. 
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a Fourier series as follows: 


w(x, y) = (y) cos — 


The Fourier expansion 


1 
- (-1) (4) 


convergent for - 3< a< 5 , is used to represent the right hand side of Equa- 


tion (2). When Equations (3) and (4) are substituted into Equation (2) and the 
Fourier coefficients on each side are equated, the following set of differen- 
tial equations are found for wp(y) for the odd values of n only:* 


For n even, the w,(y) are identically zero because of the symmetry of the 
loading and the boundary conditions. For n odd, the dimensionless variables 
u and B are defined for convenience as follows: 


man afb 


and we let 


(4)= (4) 


Then Equation (5) can be written 


* Primes indicate differentiation. 
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Now substituting the Fourier series of Equation (3) into the well-known 
formulas [1] for the stress resultants, i.e., for moment, shear and edge re- 
action distributions (see Figure 2), leads to Equations (9) through (18) below. 


M(x,y) = 


M(x,y) = - 


2 n7 (11) 
Q(x,y) = zy ain sin. (12) 


x 


ni 


Q (x57) at 


aM 
4 (13) 
Y 


(azy (2%) » (4) sin (15) 


n=1,3,. 


The concentrated corner reaction forces are found from Equation (14) 
above through the relations (cf. [1] page 92): 


R(a/2, 0) = <2 My 


(a/2, 0) (17) 


R(a/2, b) = 2 — 


(a/2, b) (18) 


The edge y = 0 is assumed simply supported while the edge y = b is assumed 
free. The boundary conditions are then (using Equations (10) and (16)): 
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The solution of Equation (7), which also satisfies the boundary conditions 
(19 a,b), can be expressed in the following form: 


£4) 
+ 


Ae 


(21) 


By substituting Equation (21) into Equations (20 a,b) there results two linear 
equations for coefficients A, and Dy): 


(22) 


(23) 
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2 sink nu + COSA 

(isu) Cosh nu 
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An and Dy, can be found explicitly by Cramer’s Rule, but a better method for 
purposes of computation is to express them as a series in ascending powers 
of e~nu, The leading terms in these series have the form 


( nw) (1-4) + [2 


(34a) 
+ O( ) 


The form of Equation (21) is superior to the conventional expression* em- 
ploying hyperbolic sines and cosines. The latter is unsatisfactory particu- 
larly for computing numerical value from the series for distributed edge 
reactions, where more than just the first few terms are included in the sum. 
For w,(y) and its derivatives, in the interval 0 < y/o S 1, the expression 
within the curley brackets on the right hand side of Equation (21) remains 
bounded as n increases. Each of the separately computed terms, which must 
be summed to form this expression, are also bounded. These observations 
can be easily verified by using Equations (24) and (25). Hence, because of the 
form we have used, it follows that for any value of n there will be no loss of 
significant figures in the sum through subtraction of large magnitudes. How- 
ever, if the hyperbolic sines and cosines are used, these quantities become 
exponentially large with increasing n and cause the number of significant 
figures obtained in the sum to decrease (even to zero!) as n increases. This 
behavior would seriously affect the automatic computation program described 
in the following section. 


Computation of Series 


a) Using Equation (21) with Equation (9) through (18), the deflections and 
stress resultants throughout the plate can be calculated. In the following 
formulas, use is made of sums in closed form for certain simple trigonomet- 
ric series; these are given in Appendix A below. The series for plate deflec- 
tions is 


Mahi 
-~nuy/b 
[ An t2 + (/- Dr )nuy/b je f cos 


(26) 


* See, for example, [1] Section 42, Equation (f). 
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The stress resultants of greatest numerical interest for engineering pur- 
poses were taken as listed in Table I below, and explicit series for these 
quantities follow in Equations (27) through (36). The uniform convergence of 
these series can be demonstrated with the use of Equations (24) and (25). 
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Ofx<a/2 


y*0,y=#b 


TABLE I 


The physical quantity is obtained by multiplying the computed number on the 
right by a dimensional factor; the dimensionless ratios were chosen to agree 
with Timoshenko ( (1), Section 42, particularly Table 28). 
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| [An + D, 2= (/-D,)(/- f 
(33a) 


[ An Dy — 


L A, +2) + (/- Zu + 


In Equations (34) and (36) above, fo(x) stands for the rapidly convergent 


q series derived in Appendix A below. 7 
ee b) The first nine* quantities listed in Table I were computed by an IBM + 
Electronic Data Processing Machine, Model 701, which was programeed ac- 4 
cording to the Speedcode I manual. This program is capable of computing 
* Concentrated corner reactions R were not computed individually as such, 
but were obtained from the computed twisting moments by use of Equations 
(33) and (34). 
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numerical values of the functions represented by the right hand sides of 
Equations (27) to (32) and (34) to (36) for arbitrary values of u,w# andk. For 
predetermined values of these three parameters, a single run of the program 


computed simultaneously the nine functions at station values cs or M at 


intervals of 0.1 in the range 0 to 1.0. 

In order to obtain uniform accuracy in all nine functions, the program 
made the following convergence test for the infinite series portion of each 
function: 


< Convergence Factor (37) 


n 


where a, is the nth term of the series, S, is the partial sum, 


(38) 


and the Convergence Factor, (abbreviated as C. F.) is a preassigned positive 
number. If the C. F. is sufficiently small, and if inequality (37) is satisfied 
for some value of n and for no smaller value, then the summation error at 
this n is found to be close to (although not necessarily bounded by) the C. F. 
The program computed successive partial sums for each series up to this 
value of n, which is generally different for different series and at different 
stations. In this way, the summation error was kept approximately uniform 
for all functions and over the entire plate. 

A study was made to determine the optimum value for the convergence 
factor. Considerations of actual machine computing time, of machine round- 
off error and of limits on certain machine subprograms (e.g. calculation of 
exponentials), as well as consideration of the desired number of significant 
figures* were factors limiting the lowest value which could be taken for the 
convergence factor. 

A series of runs for several values of u, k and u and convergence factors 
of 2x10-3, 2x10-4, 2x10-5, 2x10-6, 2x10-7, were made. An examination of 
these results showed that the important function values were computed to 
five significant figures at interior points, and to at least four significant 
figures at corner points, if the convergence factor was 2x10-5. Decreasing 
the convergence factor to 2x10-6 was seen in some cases to cause a sub- 
program error to enter; choosing a larger convergence factor reduced the 
accuracy of edge and corner reaction calculations with little reduction in 
computing time. A convergence factor of 2x10-5 was chosen therefore as 
optimum. The computing time per “run” with this convergence factor was 
eight minutes on the average. 


c) As an internal check on the accuracy of computation it appeared con- 
venient to test the condition of equilibrium of the forces acting on the plate in 
the z direction. This condition states that the sum of all distributed edge 


* The program used was of the floating decimal type, in which all numbers 
are stored in the form of eight digits and an appropriate power of 10. 
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. reactions, concentrated corner reactions and the lateral load on the plate, all 
my taken with proper algebraic sign, shall vanish; its mathematical statement is: 


g 1, a+ 0 (39) 


Over half of the plate, this equation can be written in a form suitable for 
computation as 


2x 


u 1=] 1 /* Ca, 0) ox 
| a (>) 


My (8/2, 0) Myy (2/2, d) 


(40) 


A numerical evaluation of the left hand side of Equation (40), using computed 
function values will, of course, not produce an identical zero. The actual 
computed value of the left hand side is called the equilibrium condition re- 
sidual. The ratio of this residual to the lateral load on the plate, which is 
the first term on the left hand side of Equation (40), is called the relative 
equilibrium condition residual. 

The integrals in Equation (40) must be evaluated numerically. Since the 
function values at the corners are generally computed less accurately than 
the values in the interior, the integrals from 0 to 1 are expressed by means 
of a nine point, “open end,” integration formula using the function values at 


f dx = 10-5 [uous yy - 11690 yp + 33340 y3 - 55070 
9072 


+ 67822 - 55070 + 11690 ¥g* OLS (41) 


The equilibrium condition residual and the relative equilibrium condition 
residual are tabulated at the end of each run and graphs of these quantities 
over a range of plate sizes 1/5 < b/a <5 are shown in Figure 23. Although 
the residuals are small, it is surprising to note the rapid increase of the 
residual with increasing u. This result seemingly contradicts the analytical 
result that convergence of the series should improve with increasing u. The 
explanation is seen, however, from the shape of the curves for Vx (a/2, y)/qo# 
for large u (cf. Figures 18 and 19). The slope of these curves at y = 0 ap- 
proaches closer and closer to the vertical, and the curvature likewise in- 
creases, as u or b/a increases. It is to be expected and can be shown 
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(cf. Section IV-d below), that the integration error in Formula (41), when ap- 
plied to the second integral of Equation (40), will increase markedly with in- 
creasing u. 

Despite this difficulty, the equilibrium condition test does provide an in- 
ternal check on random error in addition to other external checks on ac- 
curacy discussed elsewhere. In the experience with this program, no random 
error has ever been detected in the results of a completed run. 


Discussion of Results and Comparison with References 


a) Numerical Calculations for pw = 0.3 


The program was run with uw = 0.3 and with both k = 0 and k = 1 for some 
twenty one values of b/a in the range 1/5 < b/a <$ 5. As an aid in the design 
of structures incorporating plates, these results are presented in graphical 
form in Figures (8) through (22). The sign convention for the various quanti- 
ties described is the conventional one [1] and should be apparent from 
Figure 2. Also, the method of obtaining the results for an arbitrary trape- 
zoidal loading by a linear combination of the results for k = 0 and k = 1 should 
again be evident. 

On many of the graphs are shown for comparison dotted lines or curves 
labelled b/a =co. For k = 0, these curves correspond to values for the in- 
finite strip plate. For k = 1, these curves represent the values of y = b/3 (the 
resultant load point) for very large b/a. 

There are a number of instances where it would be desirable to know 
function values to more significant figures than the graphs could present. 
These might be to check approximate formulas which are accurate for u 
large or to check finite difference methods of solution to this or similar 
problems. A table of values for small u is therefore appended to the end of 
this paper. 


b) Poisson’s Ratio Effects 


To demonstrate the effect of Poisson’s ratio yu, the plate with length to 
breadth ratio b/a = 3/4 was chosen, both because it is a median size and be- 
cause it was chosen for previous computations |2!) (see Section (c) below). 
Some typical values of Poisson’s ratio for different materials are: 


H = 0.16 fiberglass 

= 0.20 concrete 

= 0.25 reinforced concrete 

H = 0.30 - 0.33 steel and aluminum alloys 
= 0.49 rubber 


A nearly uniform increment in the computed functions was found by selecting 
the values “ = 0, 0.16, 0.3, 0.5. The results of the computations for both 
k = 0 and k = 1 are seen in the five sets of curves of Figures (3-7) and in 
Table (5). The variation in maximum deflections, stresses* and shears over 
the range of w from .16 to .5 appears to be less than 5% of the 0.3 values. 
The deviation from the “ =_0 values is in some cases considerably greater, 
however. The treatment in [2] uses the values for uw = 0 as approximation 
for the values at near-zero “~. There is some gain in mathematical simplici- 
ty in doing this, but as our data shows, the gain seems insufficient to warrant 
the errors incurred. 
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The greatest effect of variation of Poisson’s ratio is apparent in the dis- 
tributed edge reactions and the concentrated corner reactions. The general 
result of increasing Poisson’s ratio is to induce a more uniform distribution 
of reaction pressures along the supported edge. Increasing u results in de- 
creasing the maximum magnitudes of both the distributed edge reactions (cf 
Figures 6,7) and the largest concentrated corner reaction. In Table 5, look- 
ing at R(a/2, 0), which is positive in sign because it is directed to prevent 
that corner from rising, and at R(a/2, b), which is negative in sign because 
it is directed to support the load, it is seen for wu = 0, the reaction R(a/2, 0) 
at the interior is much larger numerically than the reaction R(a/2, b) at the 
free edge. With increasing yw, the tendency of the interior corner to rise 
from its supports is reduced and so R(a/2, 0) decreases numerically. The 
concentrated reaction at the free edge at the same time assumes more of the 
load, and so R(a/2, b) increases numerically (or becomes larger negatively). 


c) Comparison of Results with References 


The numerical values given by Goriupp [2] were limited to the case b/a = 
0.75 and yw = 0; case (b) of [2] corresponds to K = 0, while k = 1 is a linear 
combination of cases (b) and (c). A comparison with the values given by the 
IBM program for these same parameter settings revealed several discrepan- 
cies, which were in all cases traced back to errors in computation in [2], 
when these errors were corrected (Goriupp’s formulas computed correctly) 
agreement with the computer results was obtained to at least the first two 
significant digits. The incorrect values are found on page 95 and in Table 6, 
page 97, of [2 ; these and the corrected values are as follows: 


Goriupp Listed Corrected 
Notation Value Value 


(1/2, A) -01690 .01239 
gy (O, 0) 19663 


(1/6, 0) 017597 06736 
Gy (1/3, 0) 611215 -03787 


p (1/2, A) -03380 .02)78 


TABLE 2 


In (1) , page 219, Table 28, and in (3) , cases (7) and (8), data is presented on 
maximum deflection and maximum moment as taken from the work of B. G. 
Galerkin. Agreement with these results are obtained here to two significant 
digits. 


d) Formulas for Large u 


The exponential series representations of A, and Dy, Equations (24) and 
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(25), can be used to obtain approximate, closed form representations for the 
Fourier expansions. In these “asymptotic” representations, the neglected 
terms vanish exponentially as u tends towards infinity. To provide an exter- 
nal check on the computation, asymptotic formulas were obtained for the con- 
centrated corner reactions R(a/2, 0) and R(a/2, b) and for the distributed 
edge reaction V, (a/2, b). 


1. The asymptotic forms of Equations (33a), (33b) are 


2 2-a)(1 16.828 


+ O(ue™") 


R(a/2,0)q, 
Writing this as 
R(a/2,b)/a, a’ a ry (1-k) 


R(a/2,0)/a, a k/b + 


the values of rj, rg, rg, rq as functions of w are given in Table 3. 


2027712 083333 135689 
.018773 2017508 -075000 122120 
010392 033922 066667 «108550 


.002519 2049341 .058333 


-.011877 2077536 20):1667 20678), 


TABLE 3 


The influence of Poisson’s ratio upon the magnitude of the concentrated cor- 
ner reactions is considerable, as the Table shows. Values of the concentrat- 
ed corner reactions at the free edge for uniform loading (k = 0) agree with 
those given by Holl [4] , pp. 22-23, particularly the curves for “corner loads” 
in Figure (5) for large b/a ratio. Similar agreement is found at the other 
corners with the values in Timoshenko ({1 , Table 5, page 133 and Table 13, 
Page 142). In fact, values for n@2 in the latter table agree with those found by 
Equation (57), for k = 1 and # = 0.3, to two significant figures for all plate 
ratios b/a greater than 1.2, that is u > 3.8. 


2. The asymptotic form of Equation (31) for the interval 0< y/b <1 is 
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1 2 
V_(a/2,y)/a,a = - = (1-k y/b)+ - 
x 2 ie N@1, 3,5 000 


enu(l-y/b) (2-1) - nu(1-y/b)] k e 
nu 


(46) 


At the corner y = b, the infinite series in Equation (58) can be summed, 
using a formula in Appendix A, to give 


V,(a/2,b)/q,a = - 5(1-K) (47) 


3% 
and this can be written as 


Some values for v, and vg as functions of w are given in Table 4. 


0.500000 0.180918 
0.420468 0.149695 
0.350000 0.122119 
0.28636) 0.097860 
0.178572 0.058152 


TABLE 4 


To check the integration error mentioned in Section III-c above, Equation (46) 
is integrated over the interval 0 < y/b < 1 to obtain the following asymptotic 
form: 
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[ V (a/2sy)dy = (1-k/2)* + (1-k) 


- 2( 14) (2-4) uel 


+ ~ + | (49) 
Qy (0,0) 
0.35309 
0.34902 


0.34575 


0.34150 


0.22830 
0.22681 


0.22562 


0.22406 


CORNER REACTIONS AND CENTRAL SHEAR _ 
FOR VARIOUS POISSON’S RATIOS; b/a = 0.75 k = 0 andk = 1 


3 t 
3 a 
T - 0.08018 0.06737 
- 0.004034 | 0.05580 
- 0.009356 | 0.04704 
4a T - 0.013629 0.03964 
TABLE 5 
~ 


1196-40 EM 2 April, 1957 


APPENDIX A 
} q SUMMATION OF ELEMENTAL SERIES 


a Upon repeated integration of Equation (3), the following Fourier series 


are obtained which are uniformly convergent in the interval - +4 x ; and 


represent there polynomials in ( =): 


2 
of (1A) 
N*1, 3,006 né 
nel 
: (-1) 2 Y 2x 
n@1,3,... n3 (3 
(-1) 2 7* 1 
q -3 @)] (3A) 
niyx 
cos —— 5 - 6\—) + 
3,000 153% G (=) (4A) 
wo These same series, with the sine replaced by the cosine, or vice versa, are 
+ { no longer Fourier expansions of polynomials. Such a series is fp(x), defined 
as: 
= 
- 
<3 a This Fourier series converges too slowly for practical computational pur- 
j poses. A rapidly converging series representation for f9(x) is found as fol- 
lows. Let = [1 = . Then Fo(s) = fo(x) and 
n@1,3,... 
Now, by [5] page 350, 
1 > = cos = log tan 
az n=1,3, eee 
1 2 7 (7A)* 
’ * cf[6], Formula #784. The Bernoulli Numbers B,, here are defined in (7) . 
B2 = 1/6, Bg = -1/30, Bg = 1/42,... 
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Then integrating Equation (7A) 


z 


there results the following rapidly converging expansion for f2(x) which is 
used in the computation of Equations (34) and (36): 


3,006 2ent2 n(2n+1)! . 


a< Bea 


Harmonic series which have simple sums are evaluated by means of the fol- 
lowing formulas from [7] pp. 236 to 240: 


= 

5,7 9A 
n-l Pp ( ) 
(-1) 


2p+1 (10A)* 


N™1, 3,006 n 


The following harmonic series, however, can be summed by use of an 
identity for the Riemann Zeta Function ¢ (2) (ct. [8] p. 269) 


oo 


2. cw (114) 


n®1,3,... n 


For 2 an even integer, this formula is identical with Equation (9A); for 2 an 


odd integer, the following six-place numerical values for the zeta function are 


useful: 


gt (3) 1.202057, £(5) = 1.036928 
(7) 1.008349. 


* E,, are the Euler numbers: EQ = -1, F4 = -5, Eg = -61,... 
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NOTATION 


Rectangular Coordinates 

Thickness of Plate 

Intensity of Continuously Distributed Lateral Load 
Loading Contour Factcr: q = qo (1 - k y/b) 
Modulus of Elasticity in Tension and Compression 
Poisson’s Ratio 

Flexural Rigidity of Plate = ahaa 

Plate deflection - 


Bending moments per unit length of sections of a plate 
perpendicular to x- and y-axes respectively 


Twisting moment per unit length of a section of plate 
perpendicular to x-axis 


Shearing Forces parallel to z-axis per unit length of 
sections of a plate perpendicular to x- and y- axes re- 
spectively 


Edge Reaction Forces parallel to z-axis per unit length 
of sections of a plate perpendicular to x- and y-axes 
respectively 


Concentrated Reaction Forces parallel to z-axis 


Side dimensions of rectangular plate 


= nuy/b = nt y/a 
A positive integer used for summation index 
A, D, Series coefficients 


wply) = Wy (B) Fourier coefficient in expansion of w(x,y) 
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PLATE BUCKLING IN THE STRAIN-HARDENING RANGE 


Geerhard Haaijer@ 
(Proc. Paper 1212) 


1. SYNOPSIS 


The application of plastic design to continuous frames constructed of 
wide-flange shapes, imposes more severe limitations on the geometry of 
these shapes than conventional elastic design. In regions where yielding 
starts first, the flanges must be able to sustain strains considerably larger 
than the yield strain without the occurrence of local (plate) buckling. 

With this practical application in mind, the problem of buckling of steel 
plates compressed beyond the yield strain is treated in the present paper. 

In the strain-hardening range the material is considered to be homogeneous. 
However, because of the yielding process the material cannot be expected to 
remain isotropic. Therefore, general expressions for the buckling strength 
are derived assuming the material to have become orthogonally anisotropic. 

Orthogonal anisotropy in the case of plane stress is expressed mathe- 
matically by stress-strain relations involving five moduli. Numerical values 
of the moduli are estimated from the incremental theory of plasticity taking 
the second invariant of the deviatoric stress tensor as the loading function. 
The influence of initial imperfections is taken into account through proper ad- 
justment of the values of the moduli. In selecting these values due considera- 
tion is given to the results of buckling tests. 

In the yielding range the average strain in the direction of loading is be- 
tween the strain at which yielding starts and the strain at the beginning of 
strain-hardening. For this case the material is considered to be partly 
elastic and partly strained up to the strain-hardening range. 

Finally, theoretical estimates are compared with test results. It is con- 
sidered that the theory adequately describes the behavior. 


Note: Discussion open until September 1, 1957. Paper 1212 is part of the copyrighted 
Journal of the Engineering Mechanics Division of the American Society of Civil 
Engineers, Vol. 83, No. EM 2, April, 1957. 


a. Fritz Eng. Lab., Dept. of Civ. Eng., Lehigh University, Bethlehem, Pa. 
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2. INTRODUCTION 


Presently used steel wide-flange shapes are proportioned such that no 
local buckling occurs within the elastic range. Consequently they can safely 
be used for structures in which the design is based upon theoretical first 
yield as the limiting condition (conventional design). However, design based 
upon ultimate strength (plastic design) imposes more severe requirements on 
the sections with regard to local buckling. The structure will reach its full 
ultimate load only if those parts where yielding starts first, can undergo suf- 
ficiently large deformations. For framed structures constructed of wide- 
flange shapes the flanges at the above mentioned locations must then be able 
to sustain strains considerably larger than the yield strain. Consequently the 
flanges should be proportioned such that local (plate) buckling does not occur 
under this condition. 

In order to solve problems of plate buckling the relationships between the 
increments of stresses and strains due to the deflection of the plate out of its 
plane must be known. Within the elastic range the assumption that the mate- 
rial is isotropic and homogeneous leads to predictions which are in good 
agreement with test results.(1) A satisfactory transition curve for the range 
from the elastic limit stress to the yield stress can easily be obtained by 
applying Bleich’s semi-rational theory to an effective stress-strain curve.(2) 

During the yielding process the material is heterogeneous. Yielding takes 
place in so-called slip bands and the strain jumps from its value at the elas- 
tic limit to that at the beginning of strain-hardening.(3) When all the materi- 
al has been strained to the strain-hardening range the material again becomes 
homogeneous. In the strain-hardening range, stress-strain relations of dif- 


ferent theories of plasticity could be applied. Such theories can be divided 
into two groups: deformation or total stress-strain relations and incremental 
stress-strain relations. 

Bijlaard(4) was first to apply deformation stress-strair. relations to the 
plate buckling problem. The theory was developed further by Dyushin(5) and 
modified by Stowell. (6) Incremental stress-strain relations were applied by 
Handelman and Prager.(7) An extensive met | of stress-strain relations in 


the plastic range has been made by Drucker.(8) Although the necessity for an 
incremental type of mathematical theory of plasticity has been shown, the 
results of plastic buckling tests on aluminum plates are well correlated by a 
deformation theory and bear no resemblance to predictions of incremental 
theory. Onat and Drucker(9) investigated the influence of initial imperfec- 
tions on torsional buckling for a simplified model of a cruciform section. 

For this case the paradox appears at its worst. Onat and Drucker showed 
that incremental plasticity leads to proper results when unavoidable initial 
imperfections are taken into account. 

All theories of plate buckling in the plastic range imply orthotropic be- 
havior of the material. This assumption seems to be very reasonable. 
Therefore, in the present report general expressions for the buckling strength 
of orthotropic plates are derived from general stress-strain relations in- 
volving five moduli (Chapter 3). Tests of steel tubes under combined com- 
pression and torsion showed that the behavior of the material is well de- 
scribed by the incremental theory with the second invariant of the deviatoric 
stress tensor as the loading function.(10) Consequently, this theory is used 
in order to obtain values of the moduli of the general stress-strain relations. 

Generalities on stress and strain and incremental stress-strain relations 
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are summarized in Chapter 4. The influence of initial imperfections is illus- 
trated in Chapter 5. From the results of coupon tests numerical values of 
the moduli are then obtained in Chapter 6. The influence of initial imperfec- 
tions is taken into account through adjustment of the values of the moduli. 
Combining the results of Chapters 3 and 6 gives numerical solutions of the 
plate buckling problem which are compared with test results in Chapter 7. 
In summary, then, the objective is to predict the strain at which buckling 

occurs in steel plate elements when the strain has exceeded the elastic limit. 


3. Buckling of Rectangular Orthotropic Plates 


3.1 General 


Consider a rectangular steel plate taking the center plane of the plate as 
the x - y coordinate plane. Compressing the plate in the x - direction into 
the strain-hardening range may affect all deformation properties of the ma- 
terial. Hence the tangent moduli, Ex and Ey in the x - and y - direction re- 
spectively, are probably different. The same may hold for the coefficients of 
dilatation, vy and vy in the x - and y - direction. The shear modulus, G, 
may also be affected. 

Thus 


normal strain 
shear strain 
normal stress 
shear stress 


Then the relations between the increments of strains and stresses can be 
written as follows: 


dé,=— 


de, 
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If equations (3.2) are valid for the entire cross-section the expressions for 
the bending and twisting moments in terms of the deflection, w, in the direc- 
tion of the z-axis become 


ow 3.3 


(3.4) 
(3.5) 


t = thickness of plate. 


The condition that the bent position is in equilibrium can be expressed by 
the following differential equation: 


D, + 2H + D, dw (3.6) 


Duy + Dyy + 4G, 


The derivation of these equations may be found in the pertinent litera- 
ture.(11) Only if H2 = DyDy, an assumption made by Bleich,(2) can solutions 
of this differential equation be easily obtained. 

If the plate is initally perfectly plane the value of go, at which bifurcation 
of equilibrium occurs (the plane and the bent position are both equilibrium 
positions) is determined by equation (3.6). The condition that both the plane 
and the bent position are equilibrium positions can also be expressed in terms 
of work. The additional work done by the external forces due to bending of 
the plate must equal the change in internal energy of the plate. 
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This yields the following integral equation 


When external restraints are provided to the plate the right-hand side of 
equation (3.7) has to be supplemented by additional terms expressing the work 
done by these restraints. 

By assuming an appropriate deflection surface, equation (3.7) gives an ap- 
proximate solution. The degree of approximation depends upon the correct- 
ness of the assumed deflection surface. In any case the result will be 
conservative. 


3.2 Plates with One Free Edge 


For a rectangular plate with the loaded edges x = 0 and x = ¢ hinged, edge 
y = 0 restrained against rotation and edge y = b free (Fig. 1) the following de- 
flection surface is assumed 


| sin (3.8) 


= | 


The ratio B/A depends upon the amount of restraint. In the case of elastic 
restraint, where X = moment per unit length required for a unit rotation 


yo (3.9) 


Deflection surface (3.8) is similar to the one used by Lundquist and 
Stowell.(12) It can be shown(20) that better results are obtained with equa- 
tion (3.8) if the following values for aj and ag are used 


For 0< 8<0.1 aj = -0.7 
ag = 0.2 

and for B= ay = -1.10 

ag = 0.54 
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Substituting w in equation (3.7) and integrating gives ; 


-(D +D Bc,+BCs 

1 46. 1+ Cr 3.10 
ar" 1/2 + 2/5 a; + 1/3 ay 


Co = 1/5 41/3 a, + 1/7 + + 1/4 aja, + 1/9 
C3 = 4+ 12a, + 144/5 a, + 12a; + l6az + 36a) a2 

Cc, 1 + 2a; + 3a9 

= 2/3 + 2a; + 1/5 + + + 12/7 af 

Ce = + a, + a>) 


4/3 + 3a, + 1/5 (9a; + 16a,) + + 16/7 


8 The minimum value, ¢¢;, of 0x is obtained for (/b given by 

= 

¥ 28 + B Cs Dy 

7 In the limiting cases when the edge y = 0 is hinged or completely fixed 

a equation (3.10) reduces to 

’ a.) Edge y = 0 hinged (8 = 0) and f = L 


For a long plate the first term can be neglected and 


7 i b.) Edge y = 0 completely fixed (8 = 
The minimum value, of 0x is obtained when the half-wave length 
satisfies 


| 
ik 
i, 4 
wen 
{ 
4 
a «x (3.12) 
= 
| 
id 
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(3.14) 


= [0.769VD.D, = 0.270 (Dw + 1.712 (815) 


3.3 Plates Supported Along All Four Edges 


The loaded edges x = 0 and x = ¢ are hinged and the edges y = + d/2 have 
equal restraint against rotation (Fig. 2). For this case the following deflec- 
tion surface is used(13 


+(A+B)- cos sin (3.16) 


The ratio B/A depends on the amount of restraint. For elastic restraints 
with X = moment per unit length required for unit rotation 


B_ ed (3.17) 
pe a. 


Substituting w from equation (3.16) in equation (3.7) and integrating gives 


4+Bc,+Be, 


2 £+(c.+2,)B+ Bes 
/ 


+ 


+ Bc, + Bx, | 
4 + Bc, + B%, | 


c, = 0.09472 
cy = 0.00921 
cg = 0.04736 
c4 = 0.01139 


The minimum value of 0, is obtained for ¢ given by 


(3.19) 
d D, +(C,+ &)BtBT, 


= ASCE HAAIJER 
P 
D, 
Then 
" 
ig 
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x 12 Fy Ls 
2 i 
1 + (Dy + Dn) 
(3.18) 
4 
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4 
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Fig. 3 - Simplified stress-strain curve 
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In the limiting cases, when the unloaded edges y = ¢ d/2 are hinged or 
completely fixed, the minimum values of ¢x are 


a.) Edges y = + d/2 hinged (8 = 0) 


(3.23) 


In the following chapters values of Dx, Dy, Dxy and G; will be determined. 
On substituting these values in the above general expressions, numerical 
solutions to the local buckling problem will be obtained. 


4. Incremental Stress-Strain Relations 


4.1 General 


It is commonly assumed that yielding occurs whenever some function of 
stress, f (94j),* equals some number, k. If the material is originally isotropic 
this yield condition is independent of the orientation of the coordinate system. 
In this instance f must be a function of the stress invariants. An extension of 
the yield function is obtained by assuming the existence of a loading function, 
f (ij), which depends upon the state of stress and strain and the history of 
loading. For ideally plastic materials plastic flow occurs whenever f equals 
some number k. For materials exhibiting strain-hardening plastic deforma- 
tions occur when the loading function exceeds k. 

Prager(15) proved that, if 


* Tensor notation is used in referring to generalized stress and strain. 
Cartesian coordinates x1, x2, and x3, corresponding to the x, y, and z axis 
of engineering notation are denoted by letter subscripts i, j, k, 1, which 
take the values 1, 2, and 3. Thus the nine components of stress and strain 
tensors are represented by single symbols 0 jj and €jj respectively. Re- 
peated subscripts indicate summation. See e.g. Ref. 14. 


ASCE 
Ag 
= (4) [ YA, + +D,, + 4G, (3.20) 
where 
/ D, (3.21) 
7 b.) Edges y = + d/2 completely fixed (B = -~) a, 7 
~ (4 +1237 (D..+D..) + 
re tay 4.52 LD, (hy (3.22) 
+ 4.943 G, 
where ; 
~ 
‘1D : 
a Vy 
| | 
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1. a loading function exists, and 
2. the relation between infinitessimals of stress and strain in linear, 


the only permissible stress-strain relation for strain-hardening material 
when loading is 


(4.1) 


and when unloading is 


(4.2) 


rf = plastic component of strain € ij and F and f are functions of 


stress and strain. The geometric proof of Prager’s stress-strain law (equa- 
tions 4.1 and 4.2) is also included by Drucker in his survey of stress-strain 
relations in the plastic range.(8) 

As no information was available concerning the actual behavior of steel in 
the strain-hardening range a few tests were carried out on combined com- 
pression and torsion of steel tubes.(10) The tubes were compressed into the 
strain-hardening range and then subjected to torsion while keeping the axial 
load constant. It was found that for this particular loading path the behavior 
is very well described by Prager’s incremental stress-strain relations tak- 
ing f = J2, where Jo is the second variant of the deviatoric stress tensor.* 

Although these tests are by no means a general verification of this theory 
they give some indication of its possible validity. In view of these results 
and on account of its simplicity, the loading function f = Jo will be applied in 
the following derivations. 


where € 


4.2 Loading Function f = J2 
Applying the loading function f = Jg to equations (4.1) and@.2) gives 


(4.3) 
when dJg > 0 
and 


(4.4) 


when dJ» <0 
The increments of the elastic components, € ij’ of the strains are given by 


Hooke’s law 


* The state of stress, with components 0jj, can be split into two parts: a 
uniform tension (or compression), s, and another state of stress, with com- 
ponents s;jj, having the same shear stress but zero mean normal stress. 
The latter is called the deviatoric stress tensor. Thus 04j = sjj + 6 ij with 
= 1/3 The second invariant is given by J2 = 1/2 84; 
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(4.5) 
modulus of elasticity 
v = Poisson’s ratio 
ij = Kronecker delta, defined as unity for i = j and zero for i # j 
For the case of plane stress (07 = T xz = Tyz = 0) the stress-strain 
relations, written in unabridged form, are 
dJ, = (20,- 0) do, + (2a,-¢,)da, + 2TdT >O (410) a 
/ 
de, = tde (4.11) 
x E x E y 
Vv 
dd, 
dez =- (do, + doy) 
dt 


dJ, <0 


The function F can be obtained from the results of a simple coupon test for 


which = T =doy=dT=0 


us? 


q 
when 
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Denoting 


(4.15) 


F is defined by equation (4.6) as 


-+] (4.16) 
V2 t 


Because of unavoidable initial imperfections the above derived stress- 
strain relations cannot be applied without modification to the local buckling 
problem. After investigating the influence of initial imperfections on two 
simplified models in the next chapter effective stress-strain relations for the 
strain hardening range of steel will be derived in Chapter 6. 


5. Influence of Initial Imperfections 


5.1 General 


A perfectly plane plate will remain plane if it is subjected to loads acting 
in its center plane which do not exceed the corresponding buckling loads. In 
the case of longitudinal loading in the x - direction producing a state of stress 
with o, as the only component, this state of stress will remain unchanged up 
to the point when buckling occurs. Consequently the buckling stress can be 
obtained from stress-strain relations (4.6) to (4.9) with Oy=T= 0. 

However, the buckling strength of actual plates with unavoidable imperfec- 
tions does not agree with the predictions for perfectly plane plates. The main 
reason for the discrepancy seems to be equation (4.9) which predicts elastic 
behavior with regard to the superimposed shear stresses. 

Applying a simplified stress-strain diagram to a simplified model of a 
cruciform section Onat and Drucker(9) have shown that small unavoidable 
imperfections may account for the difference between predicted and actual 
behavior. Apparently the influence of imperfections on sections which fail by 
torsional buckling is completely different from those which fail in bending. 
The latter case has been investigated by Wilder, Brooks and Mathauser.(16) 

In the following, this difference in behavior will be illustrated for simpli- 
fied models which buckle in the strain-hardening range. The applied simpli- 
fied stress-strain curve with n = E/Et = 40 is shown in Fig. 3. Reasons why 
the compressive stresses can exceed the yield stress, 0,, will be discussed 
in Chapter 6. 


5.2 Simplified WF Column 


The simplified WF column consists of two thin flanges of equal area sepa- 
rated by a web of infinite shear stiffness and negligible area (Fig. 4). Instead 
of a true initial imperfection, the deflections at the beginning of strain- 
hardening (¢ = J, and € = €,) are used in the computations. 

Following the same approach as Wilder, Brooks and Mathauser the deflec- 
tion curve is assumed to be 


y= Sin (5.1) 
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Cross - section 


Fig 4 - Simplified WF column 


Pp 
CYOSS - 
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At the beginning of strain-hardening 


= ¥, Sin (5.2) 


The load vs deflection curve is found by considering equilibrium of the center 
section of the column. 

For the first part of the load vs deflection curve the strain in both flanges 
increases and the relation between average stress and deflection is given by 


average stress of both flanges 
(LAY 
y/d 

= depth of section 

= value of f at beginning of strain-hardening. 


(tangent modulus stress) 


Strain reversal occurs for 


The corresponding stress, 7., is obtained from equation (5.3) by substituting 
f = fs. 
After strain-reversal has started the load vs deflection relation is given by 


0, L2(n+!) To "+! 


Figure 6 shows curves of a/05 vs f for %+/%o = 1.2 and different values of 
{,- The figure illustrates the behavior of the column for loads corresponding 
to stresses 9 = 0;. Although the deflections start to increase more rapidly 
the load continues to increase. Therefore it is safe to use the tangent modu- 
lus load, which corresponds to ¢ = 6+, as the limit of usefulness of the column. 


5.3 Simplified Cruciform Section 


In contrast to the above example the influence of initial imperfections on 
the buckling strength of a column of simplified cruciform cross-section will 
now be illustrated. The solution of the problem as given by Onat and 
Drucker(9) can be applied without modification to the simplified stress strain 
curve of Fig. 3. 

The cross-section consists of a thin shell of constant thickness h (Fig. 5). 
The column which is loaded uniformly is assumed to fail by twisting. The 
ends are considered as providing no restraint, which considerably simplifies 
the kinematics of the problem and makes the state of stress and strain the 
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Same at each cross-section. An approximate solution for small values of X 
is given as 


(5.6) 


angle of twist per unit length 


Se = G (elastic torsional buckling stress) 
b? 
Results for the strain-hardening range of steel, n' = 46 (n = 40) are shown 
in Fig. 7. Load vs twist curves are plotted for initial imperfections b6@, = 0, 
0.01° and 0.19, b@ being the angle of twist between two cross-sections a dis- 
tance b apart. In all cases the ratio of the elastic buckling stress to the yield 
stress is five. 


It is seen that very small imperfections cause a considerable reduction of 
the column strength. A perfectly straight member would reach its elastic 
buckling load, for the case considered ¢e/¢, = 5. An imperfection at ¢ = go 
and € = €,, of b@, = 0.019 reduces the maximum load to 0m/¢o = 1.4. Con- 
sequently the application of the J2 incremental theory to a perfectly plane 
plate which fails primarily by twisting cannot be expected to correctly predict 
the buckling strength of actual plates. 

Rather than attempt to solve the buckling problem of a plate with initial 
imperfections, effective stress-strain relations are determined in the next 
chapter. It will be necessary to modify the values of Dy, Dy, Dxy, and Dy x, 
and G; such that the application to the general expressions of Chapter 3 will 
result in a correct description of the behavior of actual plates. 


6. Stress-Strain Relations for the Strain-Hardening Range of Steel 


6.1 Results of Coupon Tests 


A simplified stress-strain curve obtained from a simple coupon test is 
shown in Fig. 3. It must be borne in mind that the strain represents an aver- 
age strain measured over a certain gage length. It would be entirely errone- 
ous to assume that the local strains within the plastic range from €¢ to €o are 
equal to the average strain. Yielding of mild steel occurs in small slip 
bands. (3) Slip takes place in a “jump” such that the strain across such a nar- 
row band jumps from €f to €o9. The first slip band originates at a weak point 
in the specimen, due to an inclusion, a stress concentration or other defects. 
From there on yielding will spread along the specimen. 

This consideration leads to the conclusion that there is no material within 
the specimen at a strain between the yield strain, € ¢, and the strain-hardening 
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strain, €,. Either the material is still elastic or it has reached the strain- : 
hardening range. : 
In the strain-hardening range, € > € 9, the material is again homogeneous . 


and in this range the “Jo theory of plasticity” will be applied. In the inter- 
mediate range, € ¢ << € <€o, the specimen can be considered to consist of two 
materials. 

The results of 21 compression coupon tests carried out at Fritz Engineer- 
ing Laboratory are summarized in Table 1. The coupons were cut from the 
flanges of WF shapes and from angles. For the interpretation of the results 
of coupon tests the following must be taken into account. Coupons are tested 
continuously in a hydraulic testing machine. It has become customary in 
Fritz Laboratory to test coupons with a valve opening of the machine corres- 
ponding to a strain rate of 1 micro in./in. per second in the elastic range. It 
has been shown by Huber and Beedle(17) that the ratio of the yield stress of a 
static test (where the load settles down after each increment of strain) and 
the yield stress of a continuous coupon test is approximately 0.925. Conse- 
quently a value of the yield stress, @o, of 0.925 x 39.2 = 36 ksi will be used in 


the following derivations. 
i Stress-strain curves for the strain-hardening range as obtained from 5 Py ‘a 
. selected coupon tests have been replotted in Fig. 8. Coupons 9 and 18 repre- ioe 


ry sent the extreme cases while 5, 15, and 17 represent tésts with “average” - 
results. 
BS The average stress-strain curve for the strain-hardening range can be 
expressed by the three parameters introduced by Ramberg and Osgood.(18) 
= 14x 10-3 in./in. 
: 900 ksi* 
j 21 
2 


Equation (6.1) is also plotted in Fig. 8. 
The information now available is sufficient to determine F (J2) defined by 
equation (4.16). From equations (4.16) and (6.1) it follows that 


F = [67.357 591.6674 | 10° (6.2) 


for Jg > 1/2.0° = 432 kips?/in.4 


* The values of Eo in Table 1 are taken from Fritz Laboratory reports in 
which they are usually not given as the slope of the stress-strain curve at 
the initiation of strain-hardening but as the slope at a strain somewhat 
larger than €o. Consequently Eo as used in equation (6.1) is larger than 
the average value given in Table 1. 
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TABLE 1 
q RESULTS OF COMPRESSION COUPON TESTS 
Fritz Lab. 
Number 
l 220A-UF3 All WF section 
2 220A-UF4 oS -6 coupons taken 
3 220A-LF1 5 -6 from flanges 
4 220A-LF2 .0 5 
5 220A-LF3 
6 220A-LF4 .0 .0 
7 220A-A 8 .0 
8 220A-B 8 5 
220A-E 
220A-F 3 -0) 
220A-G +2 -5) 
220A-B2F3 
220A-B2F6 
220A-B2F7 -0 8 
205E-C14 
205E-C15 .0 
205E-C2 4 
205E-C9 
205E-C12 
205E-C13 


* Numbers in parentheses not used for determining average value 
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6.2 The Tangent-Modulus in Shear 


Consider the case where shear stresses, 7, are superimposed on a con- 
stant normal stress, 0 ,, taking oy = 0 and do, =doy=0. The relations be- 
tween the increments of stress and strain given by equations (4.6) and (4.9) 
reduce to 


de, = $F o,tdt 


(6.3) 


dz dt+4Ft7dt 


Integrating equation (6.4) gives the relationship between T and ¥ as shown 
in Fig. 9 for 0, = 36 ksi and cx = 48 ksi. The corresponding slope 


2(\+V7) 
(6.4) 


Tt 


(6.5) 


is plotted in Fig. 10. 
It is seen from Fig. 10 that the value of G; drops rapidly for small values 
of ¥. However, in the region 2,000 ksi < G; < 3,000 ksi the decrease becomes 

slower. Consequently any value in this region could be selected as a useful 
value of G; for the strain-hardening range of steel. On account of the results 
of torsional buckling tests on angle specimens presented in Chapter 7, the 
value G; = 2,400 ksi is selected as being applicable to the strain-hardening 
range. From Fig. 9 it follows that the influence of the magnitude of the nor- 
mal stress can be neglected for that part of the strain-hardening range under 
consideration. 


6.3 Bi-Axial Normal Stresses 


For regions of a plate in which cross bending is of importance the shear 
stresses are zero or very small, e.g. the center of plates supported along all 
four edges or the fixed edge of a clamped outstanding flange. 

In this case equations (4.6) and (4.7) reduce to 


[ 
LE 


de, =| F (20x- oy) doy + 


- $F (26, - 


doy 


(6.6) 


F (20,- | do, + 


Comparing equations (6.6) and (6.7) with equations (3.2) gives 


(6.7) 


F (20x-9y) (2ay do, 
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Fig: 10 The Tangent Shear Modulus 
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4 Ey = F (6.9) 
-5 F (26,-0,) (26, (6.10) 
2 
= +5 F (2ax-0,) 
to 
For a perfectly plane plate (Vy = 0) equations (6.8) to (6.11) reduce to 
(6.12) 
4 (6.13) 
| 4 
(6.14) 
(6.15) 
3 Equations (6.12) to (6.15) have been applied with different notation by 
Handelman and Prager.(7) 
“ Equations (6.9) to (6.11) are valid only if 4 
dJ, > 0 (6.16) 
or rewritten 
2do, - da, > 0 (6.17) 
2 4 and with equations (6.6) and (6.7) 
(2-v) de, (I-2v) de, > O (6.18) 
; Figure 11 shows the assumed linear strain distribution due to curvatures 
y dKx and dKy in the x - and y - directions. 
de, = dé, +zdK, (6.19) 
a é, = de, +2dK, (6.20) 
Ms where d€, and d€9 are strain increments of the central plane in the x - and 
P y - direction and z is the distance to the central plane. % : 
is The condition that all of the section is deformed plastically is obtained by aa 
1 3 substituting equations (6.19) and (6.20) in equation (6.18) a 
He 
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(2-v)(de,+zdk,) - (1- 2v)(de,+zdK,) > O (6.21) 
for t/2 


The increase of the force per unit width in the x - direction, Nx, is found 
by rearranging equation (3.2) and integrating over the thickness of the plate 


D,t (de,+ Vv, de,) 


dNx 


" 


The increase of the force per unit width in the y - direction is 


dN, = do,)dz = 


(6.23) 


= Dt (de, +V, de,) 


However, no external forces are applied in the y - direction, thus 


dN, =0 


de, => dé, 
Substituting equation (6.25) in (6.22) gives 


dN, = E,t de, 


The plasticity condition, equation (6.21), then becomes 


[2-v+ V, |de, +(2-v) zdk, -(I-2v)adky > 0 (6.27) 


. 
for - t/24 z4 t/2 
{ If the neutral zone between loading and unloading zones is at z = t/2, 1 
equation (6.27) gives 
: 
4 
2-V+ (1- 2V) Vy 
Obviously de > 0 only if 
uy 
T 4 dk, < dk, (6.29) 


Checking the plasticity condition (6.27) for de 1 given by equation (6.21) shows 
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that condition (6.27) is not violated if (6.29) is satisfied. 
If the neutral zone is at z = - t/2 equation (6.27) gives 


da t/2l(i-2v) - (2-'v) dkx} (6.30) 


2-V+(1-2V) Vx 


and d€, > 0 only if 
dk, > dk, (6.31) 


The plasticity condition (6.27) is not violated if equation (6.31) is satisfied. 
From equations (6.28) and (6.30) it is seen that 


dé,=0 


and consequently according to equation (6.26) 


Furthermore dJ9 = 0 for the entire cross-section. Thus 


do, = + do, (6.33) 


for an initially plane plate with ay = 0. 

In this case, since bending is not accompanied by an increase in axial load 
the influence of initial imperfections will be the greatest. Suppose biaxial 
loading starts at ox = 0x, dy =0, €x = €x, €y=€y- Then it follows from 


equation (4.10) with equation (6.33) that initially 


dJ; = (6.34) 


If during biaxial loading the ratio of do, and doy is taken according to equa- 
tion (6.33), then integrating equation (6.34) gives 


L* a? (6.35) 


Applying equations (6.8) to (6.11) to the computation of the moduli D,, Dy and 
Dxy as defined by equation (3.6) gives: 


(6.36) 


(6.37) 


= 
q 
for 
dk, = 45¥ dk, (6.32) 
av 
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Dey = Dx = 85-7? (6.38) 
where 
R= (6.39) 
$F of (6.41) 


The results are shown graphically in Fig. 12 forgo : = 36 ksi and o - = 54 


ksi. From these results it can be concluded that the influence of the magni- 
tude of o > is negligibly small. It is not obvious, however, which values of D,, 


Dy and Dy, should be selected as being applicable in the strain- hardening 
range. Fortunately, compression tests on wide-flange sections presented in 
Chapter 7 reveal that the ratios of the half-wave length of the buckled shape 
over the depth of the section for the cases of web buckling are 0.55 and 0.54 
(Tests D4 and D6 respectively). According to equation (3.21) this ratio is 


equal to 4/ Dx/Dy if the small restraining effects of the flanges are neg- 
lected. Figure 13 shows the influence of oy/a* on 4/Dx/Dy - On account 


of the results of the web buckling tests the values of Dx, Dy and Dxy corres- 
ponding to oy/o z= 0.34. Have been selected as applicable to the strain- 
hardening range. 

Thus 


D, =  3,000ksi 
32,800 ksi 
Dxy = Dyx = 8,100 ksi 


" 


7. Comparison with Test Results 


For the selection of applicable values of Dy, Dy, Dxy, Dyx and Gt use has 
been made of the results of local (plate) buckling tests. The results are sum- 
marized here and are presented in more detail in another paper. (19) 


7.1 Compression Tests on Angles 


A number of compression tests on angles were performed with the purpose 
of checking the theoretical estimates developed above. Angle specimens have 
better known boundary conditions than WF sections and therefore give a more 
positive check. When buckling torsionally under the action of an axial load, 
the flanges of the angle act as two plates each with one free and one hinged 
edge, the heel forming the hinged edge. The loaded ends of the column were 
fixed against rotation in the testing machine. The dimensions of all specimens 
are given in Table 2. Besides the longitudinal strains at the flange tips and 
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TABLE 2 


DIMENSIONS OF ANGLE SPECIMENS 


Thickness 
t (in.) Material 


Annealed 


As-Delivered 


3 
4 2L (in.) | b (in.) 
A-21 25.0 4.87 0.383 12.70 | 5.14 
- A-22 25.0 4.79 0.381 12.60 | 5.21 Be. 
A-31 17.9 3.27 0.370 8.85 15.48 | a! 
a A-32 17.9 3.28 0.374 8.79 | 5.46 re 
q 
4 A-41 12.5 2.31 0.377 6.13 | 5.41 ie 
A A-42 12.5 2.34 0.371 6.36 | 5.35 oS 
4 
A-33 17.5 3.30 0.378 8.73 | 5.30 
A-51 21.2 4.07 0.380 | 10.70 | 5.21 
oa 
¥ 
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the heel, the rotation of the center section was measured. From the rotation 
measurements the critical average strain was determined. The results of 
the angle tests are summarized in Table 3. 

When all of the material is strained into the strain-hardening range 
(€ cr > €o) the theoretical solution is given by equation (3.12), the length of 
the angle specimen being 2L. A solution for the yielding range (€ f< € cr<€o) 
can be derived if the following three assumptions are made: 


1. The material is elastic up to the yield stress. 

2. In the yield range the specimen is partly elastic and partly strain- 
hardened. 

3. The strain-hardening zones initiate at both ends and move toward the 
middle. 


From the first assumption it follows that, if 0 ¢, < %o the critical stress is 
given by 


When 6c, obtained from equation (7.1) exceeds 9 yielding will have com- 
menced. From the second assumption it then follows that the middle section, 
being still elastic, is practically rigid compared with the yielded zones. As- 
suming that only the latter will deform results in the following expression for 
the buckling stress 


where 
£ L = length of each yielded zone 


The corresponding critical strain is 
= + SE (7.3) 


Substituting in equation (7.2) the values of Dy and G; this equation determines 
the relationship between b/t and. For L/b = 2.65, the average value of the 
tested sections, the b/t vs £ curve has been plotted in Fig. 14 as a solid line. 
As elastic deformations have been neglected in equation (7.2) b/t =o for 
€=0. For this case b/t = 20.7 which is found from equation (7.1) by taking 
¢ = 36 ksi and L/b = 2.65. Knowing the rigid plastic solution and the point for 
£ = 0 of the elastic-plastic solution the latter has been sketched in Fig. 14 as 
a dotted line. The elastic-plastic solution of b/t vs £ with equation (7.3) 
gives €cr as a function of b/t for the range € < < The complete 
theoretical curves are shown in Fig. 15 for L/b =o and L/b = 2.65, and are 
compared with the test results. The theoretical curves give a good descrip- 
tion of the buckling strength. 


7.2 Tests on Wide-Flange Shapes 
In order to investigate the actual behavior of WF shapes with regard to 
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TABLE 3 


RESULTS OF ANGLE TESTS 


Type of 
Buckling 


torsional 


torsional 


torsional 


bending 


bending 


torsional 


torsional 


(1212-32 EM 2 — 
7 
3 
| 
ksi k si 
A-22 | -- 3.0 32.2| 
a A-31 | 34.9 16.5 | 
7 ag 
A-32 | 34.6 16.5 35.6 
A-41 | 35.3 -- | 
A-42 | 34.1 -- -- 
A-33 | 41.3 6.0 | 
i 
A-51 | 41.0 6.0 41.2) 
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local buckling, six different shapes were each tested under two loading 
conditions: 


a) Axial compression (Test D1, D2, D3, D4, D5, D6) 
b) Pure Bending (Test B1, B2, B3, B4, B5, B6) 


The dimensions of all WF specimens are given in Table 4 and the test results @ 
are summarized in Table 5. 

For the cases where flange buckling was predominant the critical strains 
of the flanges ws the b/t ratios and the theoretical curves are plotted in 
Fig. 16. The theoretical solution is given by equation (3.10). The results of §& 
tests D4 and D6 are omitted because web buckling occurred first and obvious- § 
ly caused premature flange buckling. Furthermore specimen B4 did not 
develop a major flange buckle but failed by lateral buckling. Therefore this 
result has also been eliminated from Fig. 16. From this figure it can be con- 
cluded that, if premature web buckling is prevented, the webs of the tested 
sections provide some restraint to the flanges corresponding to a value of 8 
from 0 to about 0.05. Comparing the theoretical values of the half-wave 
length over flange width ratio given by equation (3.11) and plotted in Fig. 17, 
with the measured values of the bending tests given in Table 5 shows that the 
theory gives a good description of the actual behavior. For the axial loading 
tests the half-wave length was obviously influences by web buckling, except 
for specimen D5 which had the smallest d/t ratio. 

For the cases where web buckling occurred first (Test D4 and D6) or 
simultaneous with flange buckling (Test D2) the critical strains are plotted vs 7m 
the d/t ratios as in Fig. 18. Comparison is made with the theoretical solu- 
tions given by equations (3.20) and (3.22) for the cases of zero and full re- 
straint respectively. Again favorable agreement is obtained. The values of 
(/d given by equation (3.21) and plotted in Fig. 19 necessarily agree with the 
experimental values because Dy and Dy were selected in view of these test 
results. 


8. Summary 


The results of the investigation presented in this paper can be divided into 
two parts: firstly the derivation of stress-strain relations for the strain- 
hardening range of structural steel and secondly their application to the plate 
buckling problem. 

Effective stress-strain relations were determined describing the ortho- 
tropic behavior of steel after it has been compressed into the strain- 
hardening range. The following values of the moduli were found to be 
applicable: 


3,000 ksi 
32,800 ksi 

Dyx = 8,100 ksi 
2,400 ksi 


" 


" 


Q 
" 


It is considered that the agreement between theory and test results (Figs. 
15, 16, and 18) justifies this approach to the problem. 

A direct practical application of the findings presented in this paper is the 
prevention of local buckling of outstanding flanges in continuous frames, in 
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TABLE 4 


DIMENSIONS OF WF SPECIMENS 


LOWF33 
8WF24 
1LOWF39 
12WF50 
8WF35 


10WF21 


width of flange 

thickness of flange 

distance between center planes of flanges 
thickness of web 

length of compression specimen 


length of part of bending specimen subjected to pure bending 


ASCE HAALJER a 
2b te d ty |L ful 
Spec. |Shape ion in in in| in | >/t¢| 
Bi Dl 7.95 | 0.429 | 9.37 |0.294 | 32] 32] 9.2 | 31.9 
gery B2 D2 6.55 | 0.383 | 7.63 |0.236 | 26| 26] 8.6 | 32.3 - 
B3 D3 8.02 | 0.512 | 9.37 | 0.328] 32| 32] 7.8 | 28.6 
; i B4 D4 8.18 | 0.620 | 11.57 | 0.351 | 32] 32] 6.6 | 33.0 . 
BS DS 8.08 | 0.476] 7.65|0.308| 32| 32| 8.5 | 24.8 
B6 D6 5.77] 0.318] 9.56] 0.232] 23] 26] 9.1 | 40.9 
2b = 
3 
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TABLE 5 


RESULTS OF WF TESTS 


Type of Buckling 


flange 


flange & web 
flange 
web 
flange 
web 
flange 
flange & lateral 
flange & lateral 
lateral 
flange & lateral 


flange & lateral 


Flange [wed [web | | 1/4 
Dl | 34.4] 8.5 | 8.5] 34.2 | 34.2 | 1.8 | 0.56 
p2 | 34.0] 13.5 [12.7] 34.0 | 34.0] 1.5 | 0.50 
p3 | 35.2 | 19.0 |19.0 | 39.0 | 39.0 | 1.5 0.46 
| 35.0 18.5 | 5.0] 36.8 | 35.4 ]1.5 | 0.55 
a ps | 36.6 | 17.0 |17.0] 38.0 | 38.0 | 2.1 | 0.56 aur 
a pe | 38.0 | 4.3 | 1.6] 33.8 | 37.2] -- | 0.54 
¥ Bl | -- 7.0 | -- -- -- | 2.4 -- ne if 
7 x B2 | -- | 23.0 | -- -- -- | 2.0 -- ans 
B3 -- 22.5 -- “ -- 2.2 -- 
BA | -- | 29.0 | -- -- -- we 
i q 
4 | -- -- | 2.0 -- 
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which the design is based upon ultimate strength. From the required rotation 
capacity of the plastic hinges the strains of the flanges can be determined. 
Fig. 16 then gives the required b/t ratio, for B = 0.01. 

These stress-strain relations could equally well be applied to other prob- 
lems involving the occurrence of shear and biaxial stresses and strains in the 
strain-hardening range of steel. 
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11. Nomenclature 


Tensor Notation 


F = function defined by equation (4.1) 


f = yield and loading function 


i,j,k, f, m - are letter subscripts taking the values 1, 2, and 3 


second invariant of deviatoric stress tensor 


constant 


mean normal stress 


components of deviatoric stress tensor 


coordinate axis 


Kronecker delta 
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components of strain tensor 


elastic strain component 


plastic strain component 


components of stress tensor 


Engineering Notation 
A,@,,42 - are constants 
B = constant 


b = width of plate with one free edge 


Cj, Cg, C3, C4, C5, Cg, C7 + - are constants 


Ex 
1 - VxVy 


Ey 
l-vp 
VyDx 
VxDy 
width of plate supported along all four edges 
modulus of elasticity 
tangent modulus 
tangent modulus in x - direction 
tangent modulus in y - direction 
strain-hardening modulus 
ratio of y over depth of simplified column 
value of f at initiation of strain-hardening 
value of f at which strain reversal takes place 
modulus of elasticity in shear 
tangent modulus in shear 
function defined by equation (3.6) 
thickness of sheet forming simplified cruciform section 
moment of inertia per unit width of plate 
constant 
curvature of plate in x - direction 
curvature of plate in y - direction 
twist of plate 
half-wave length of buckled shape 


bending moment per unit width of plate in x - direction 
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bending moment per unit width of plate in y - direction 
torsional moment per unit width of plate in 
axial force per unit width of plate in x - direction 
axial force per unit width of plate in y - direction 
exponent in equation (6.1) 
ratio of modulus of elasticity over tangent modulus 
ratio defined by equation (5.6) 
function defined by equation (6.39) 
function defined by equation (6.40) 
function defined by equation (6.41) 
thickness of plate 
deflection of plate 
coordinate axis 
coordinate axis 
deflection of simplified column 
maximum value of y 


coordinate axis 


B/A = coefficient of restraint 

angular strain in xy plane 

critical strain corresponding to ocr 

yield strain 

strain at initiation of strain-hardening 
normal strain in x - direction 

value of €y at which biaxial loading starts 


normal strain in y - direction 


value of € y at which biaxial loading starts 


coefficient determining length of yielded zone 

angle of twist per unit length 

Poisson’s ratio 

coefficient of dilatation for stress increment in x - direction 
coefficient of dilatation for stress increment in y - direction 
critical (buckling) stress 

elastic buckling stress 


yield stress 
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value of o at which strain reversal takes place 
normal stress in x - direction 
value of oy at which biaxial loading starts 
normal stress in y - direction 
shear stress 
function defined by equation (6.10) 


edge moment per unit length to produce unit rotation of edge 
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INELASTIC RESPONSE OF COLUMNS TO DYNAMIC LOADINGS 
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SYNOPSIS 


An analytical investigation of the response of columns to single compres- 
sive loadings of arbitrary time history such that the material of the column 
is strained into the inelastic range is presented in this paper. The column is 
assumed to have pinned-ends and to have a slight initial curvature or to have 
the loading applied slightly eccentric to the axis of the column. 

The modal superposition approach is used to obtain the equations for the 
elastic response while an idealized, one-degree of freedom, rigid-plastic 
system is used to obtain the equations for the ensuing inelastic response, af- 
ter the strain at some point in the column exceeds the yield strain for the 
material. Numerical integration is used to obtain solutions for both the 
elastic response and the inelastic response. A numerical example is worked 
out and the results compared to the experimentally determined response of 
an actual column. 


INTRODUCTION 


If a compressive force is applied slowly to a column, the classical stabili- 
ty analysis may be used to obtain the static stability load in the case of a 
concentrically loaded member. However, if the load is applied rapidly, the 
dynamic equilibrium of the column must be considered. As the load is ap- 
plied and the column deflects laterally, inertia forces tend to resist the 
lateral deflection and have the effect of stiffening the column, so that, if the 
load is removed before the column attains a large lateral deflection or 
velocity, the column may be able to support loads much in excess of the 
static loads for this short period. 


Note: Discussion open until September 1, 1957. Paper 1213 is part of the copyrighted 
Journal of the Engineering Mechanics Division of the American Society of Civil 
Engineers, Vol. 83, No. EM 2, April, 1957. 

* Asst. Prof., Dept. of Civ. Eng., Univ. af Washington, Seattle, Wash, 

** Associate Prof., Div. of Civ. Eng., Univ. of California, Berkeley, Calif. 
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If the column is a very long, slender member, the stresses due to a 
dynamic loading will generally remain below the yield point for the material 
throughout the entire response. This problem of elastic buckling of columns 
subjected to dynamic loadings has been treated by several previous investi- 
gators; Koning and Taub(1) followed by others(2,3,4) have investigated the 
elastic response due to a given applied loading and Hoff and others(5,6) have 
investigated the response due to a given time varying displacement of the 
end of the column. 

However, most practical columns have a slenderness ratio such that they 
buckle only after the stresses in the column exceed the yield point. Two 
previous investigators(7,8) have discussed the inelastic response in the case 
of a given time varying displacement of the end of the column. The purpose 
of this paper is to present a method of analysis for the inelastic response of 
a column to an arbitrary dynamic loading. The numerical method of integra- 
tion presented may also be used for obtaining the complete elastic response 
to an arbitrary dynamic loading in the case when the column does not yield. 
In addition the method may be adapted readily to obtain the elastic or in- 
elastic response of certain types of beamns subjected to dynamic loadings or 
to specified initial displacement and velocity distributions. 


Equations for Elastic Response 


The initial response to a dynamic loading will be elastic in all cases ex- 
cept the idealized case of a suddenly applied load equal to the dynamic yield 
strength of the column. Therefore, the first part of the analysis is devoted 
to obtaining the elastic response from which can be determined the instant 
of yielding and the initial conditions for the inelastic response. 

Consider the column shown in Fig. 1, in which y(x,t) is the deflection of 
the column at time t and y,(x) is the initial deflection of the column due to 
either initial curvature or end eccentricity, e, or a combination of these. 
The differential equation of equilibrium of this column for constant stiffness, 
EI, and mass per unit length, m, and neglecting rotatory inertia and shear 
effects is: 


Assuming, for pinned-end columns, 


y(x,t) = Un sin 


and 


YolX)= Uon Sin (3) 
=/ 


where 
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Uon Yo (x) sin dx (4) 
(7) 
and substituting into Eq. 1 yields: 


But since the terms of this series are all linearly independent, in order for 
the sum to be zero the coefficients of each term must be zero, or 


Rewriting this differential equation for the “generalized coordinates,” u,(t) 
gives, (using dots to represent differentiation with respect to time) 


=1,2,---00. 


Eq. 7 can be simplified by introducing 


2rr2 
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pit) (9) 


T= tw (i) 
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Eq. 7 then becomes 


t -An(7)] = 4 5 2, --- (12) 


The initial conditions are: 


44,2) = Yo(x) and 240.0) =O (13) 


or 


and (4) 


Eq. 12 with the initial conditions (Eq. 14) are an infinite system of linear, 
second-order differential equations for the quantities v,(t) from which the 
deflection can be obtained using Eqs. 9 and 12. The Solution of these equa- 
tions for constant A ,(t) (i-e., constant P(t) are well known (three solutions 
exist, one each for ,<1,A, = 1, andAy>1). It should be noted that the 
series of Eq. 2 will not converge to the correct moments at the end in the 
eccentrically loaded case and, in general, will converge only slowly near the 
ends in this case. Pian & Siddall(3) have given a form of the solution with 
more rapid convergence which corresponds to subtracting the Fourier sine 
expansion for a parabola and adding the equation of a parabola to this series. 


Equations for End of Elastic Response 


The equations derived above are valid until the strain in some fiber of the 
column exceeds the yield strain. The strain in an outside fiber of the column 
is given by (c is its distance from the neutral axis) 


where y(x,t) and yo(x) are given by Eqs. 2 and 3. For symmetrical initial 
shapes subjected to loads of not too short duration (>0.3Ty, where Ty is the 
natural period of vibration), the onset of yielding may be defined with suf- 
ficient accuracy by using only the 1st and 3rd modes of the response, and the 
maximum strain will be at the center. However, very rapidly applied eccen- 
tric loads will cause yielding at the end of the column if the strain at the end, 
given by the expression 


exceeds the yield point. In the case of unsymmetrical initial shapes, the Ist 
and 2nd modes of the response would generally be required to determine 
satisfactorily the location and instant of first yielding. The strain at which 
yielding occurs would, of course, depend on the rate of straining of the mate- 
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rial but within a certain range of loading, a satisfactory single value probably 
could be established for all cases. It should be noted that even though only 

one mode may be required to define the deflections, more than one mode may 
be required to determine adequately the instant of yielding. 


Equations for Inelastic Response 


As soon as yielding commences in the column, a rigid-plastic type of ac- 
tion is assumed in which a hinge develops at the point of yielding and the seg- 
ments of the column on either side of the hinge are “frozen” in the deflected 
shape which existed at the instant of yielding. Such a system is shown in 


: Fig. 2 for a hinge located a distance x, from the left end of the column. The i 
dynamic moment equilibrium condition of the left segment of the column pro- 
vides the equation 
ee = 
or simplifying 
where My is the internal plastic moment resistance of the hinge. For the 
ordinary symmetrical case of a hinge at the center of the column this 
equation becomes 


= PU ult)- My (9) 


The initial conditions for the inelastic response are obtained from continuity 
of displacements at the position of the hinge, and from continuity of momen- 
tum from the elastic motion to the rigid-plastic motion. For the case of a 

hinge at the center, these initial conditions for the inelastic response are 


(20) 


Uo inelastic = U(t;) elastic 


and 


A. 
Uo inelastic * “7° (7 ij) elastic @/) 


where t; is the time at the instant of yielding. 

Eqs. 18 or 19 for the inelastic response require an expression for the in- 
elastic moment My in terms of t. This inelastic moment for a column will be 
a function of the shape of the cross-section, the properties of the material, 
the axial force, and the curvature, which is in turn a function of the deflec- 
tion. In general, the true relationship between moment (M), axial force (P) 
and curvature (@) in the inelastic range can be approximated by one or more 
linear relationships such as 
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My = Ma(1- ) (22) 


where Ma,¢@a, and Pg are constants to be determined for the particular 
cross-section and material properties. It may also be possible to find an ap- 
proximate linear relationship between ¢ and u in the inelastic range of the 

form 


Gh) =kult)+C (23) 


which would correspond to a yield zone of fixed length. It should be noted 
here that although it was possible to put the equations for the elastic response 
(Eq. 7) in the non-dimensional form of Eq. 12, the equations for the end of the 
elastic response and for the inelastic response do not depend on these same 
non-dimensional parameters, hence the total response cannot be put in non- 
dimensional form. For this reason, the total response will be a function of 
the actual dimensions as well as of the initial deflections of the column and 


the loading. 


Numerical Integration of the Equations 


In the case when P(t) is a constant, Eq. 12 may be solved explicitly for the 
elastic deflection. A solution in terms of Mathieu functions is available for 
the case of a sinusoidal loading. It is not likely that a solution of Eq. 18 for 
the inelastic response could be found except in the special case of constant 
load and for a particular form of internal moment, My. For arbitrary P(t) 
and for more general forms of M], numerical integration of the equations 
must be used to obtain solutions of the equations. For dynamic equations 
similar to these, other investigators have used various step-by-step, colloca- 
tion and finite differences methods.(3,7,9) The method chosen here is a step- 
by-step iterative method which has been applied to dynamics equations by 
Newmark.(10) This method has the advantages that relatively large time in- 
tervals may be used, that the initial conditions enter in a simple manner, and 
that quantities such as accelerations and velocities do not lose their physical 
significance. 

The method of numerical integration used here consists of assuming that 
the second derivative of the displacement varies linearly over a small time 
interval, 4t, from the known initial value, up, to an assumed final value 
(ty + At). Integration of this twice yields the displacement at the end of the 
time interval. This value of the displacement at the end of the time interval 
can then be used in the equations for the response of the column to obtain a 
value for u(tp + At) which should check the assumed value. The assumed and 
computed values of l(t, + At) at the end of the interval are made to agree by 
iteration at each time interval. The method is portrayed in Fig. 3. The 
equations for u and u for an assumed value of ii(t, + At) are 


alt, +At)= A + ip (24) 
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and 


Ulty +7) = At) inp)? (as)? +Up At (25) 


Since either Eq. 7 for the elastic response or Eq. 19 for the inelastic re- 
sponse gives u(t) as a function of u(t) and P(t), the appropriate equation can 
then be used to compute ii(tp + At) from Eq. 25 and P(tp + At) from the given 
load function P(t). If this computed value is not the same as the value as- 
sumed for computing ii(ty + At) from Eq. 25, this new value can be used as 
an assumed value for repeating the process. If the time interval At is small 
enough the process will converge quickly to the correct acceleration. New- 
mark indicates that satisfactory convergence will be obtained if the time in- 
terval of integration is less than 1/10th the period of vibration for simple 
vibration problems. The time interval used in this investigation was less 
than 1/10th the period of vibration of the unloaded column. The convergence 
was checked using an interval half as large and was found to give the same 
results (within one percent). In some cases the maximum usable time inter- 
val may be the largest value that would satisfactorily approximate either the 
load curve or the instant at which the column yielded. 


Numerical Example and Comparison with Test Results 


Test results were available for about 80 dynamic loading tests performed 
by the authors on columns in which the response was in the inelastic range. 
The analysis presented in this paper is now applied to one of these tests. 
The column considered was an annealed mild steel tube, 1 inch square 
outside with a wall thickness of 0.11 inch and a length of 44.5 inches. The 
L/r ratio was 124, the Euler buckling load was 7.2 kips, and the fundamental 
period of free lateral vibration was 0.0198 sec. The static yield load for the 
section was 12.8 kips (32.0 ksi). The test results indicated that the upper 
dynamic yield point loading was about 23.6 kips (62.0 ksi) giving a dynamic 
yield strain of 2100 microinches/inch. The lower dynamic yield point was 
estimated to be 17.7 kips (46.5 ksi) at a strain of 7350 microin./in. Stress- 
strain properties of the material are depicted in Fig. 4. Using these dynamic 
stress-strain properties and the dimensions of the cross-section gives the 
dynamic moment-curvature-axial force relationship shown in Fig. 5. The 
relationship used for M] in the analysis is shown by the dotted lines, and is 
of the form of Eq. 22, with the constants defined as follows. From the instant 


of yielding (which is assumed to occur at $= 4.90 x 1075(1 - aa? until the 


curvature reaches the value ¢ = 11 x 107*in /in/in, the constants are taken 
as P, = 23.6 kips, Ma = 8 kip-inches, and @, = 44 x 10-3in/in/in. For 
curvatures greater than ¢ = 11 x 10-3in/in/in, they were taken as Pa = 17.7 
kips, Mg = 6 kip-inches and ¢g =o. The relationship between deflection and 
curvature for use in Eq. 22 was found to be given by an equation of the form 
of Eq. 23 as follows 


+ (u -u;) (2) 
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assuming a hinge to develop at the center, where $j and u; are the values of 
¢ and u at the instant of yielding. This corresponds to a uniform yielded 
length of 10 inches. 

The column was concentrically loaded and, since its initial shape was 
primarily a first mode deflection, only one elastic mode was used for the 
analysis. The initial first mode deflection was 0.043 inches. The loading 
and the measured and computed center deflections plotted as functions of 
time, are shown in Fig. 6. The time interval of integration was one-half 
millisecond (0.0005 second). On this figure is indicated the instant at which 
the column was found to yield in the analysis, according to the criterion 
established above. Values of the Euler buckling load, the static yield load, 
and the fundamental period of free vibration are also shown. 


CONCLUSIONS 


The method of analysis presented gave a very good approximation to the 
actual response of the column considered, even though only one mode was 
used in calculating the elastic response. This column had a dynamic loading 
applied during an interval equal to 0.6 times the fundamental period of vibra- 
tion of the unloaded column (Ty). 

The method was found to give good results when applied to any of the 
columns tested for which the duration of the loading was from 0.3 to 0.9 
times the fundamental period of vibration. However, if the period of loading 
were quite short (i.e., < 0.3 T)y) it would be expected that more modes would 
be required to determine satisfactorily the position and instant of yielding, 
even for columns containing relatively small amounts of higher modes in the 
initial shape. This is because the short duration loading would be more 
effective, relatively, in exciting the short period modes of vibration. On the 
other hand, if the period of loading were relatively long (i.e., > 0.9 Ty), it 
would be expected that a better approximation of the inelastic moment resist- 
ance and the dynamic yield point would be required to give a satisfactory 
approximation of the actual deflection, because the time during which dynamic 
yielding took place would occupy a more significant portion of the total re- 
sponse. Both of these expectations were substantiated by the results of the 
test program. If proper consideration is given to these factors, the method 
presented may be applied to columns subjected to any type of dynamic load- 
ing, subject only to the limitation that the columns considered form a single, 
reasonably concentrated, plastic hinge during yielding so that the inelastic 
response may be approximated by a rigid-plastic column. 
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Fig. Z.-- Inelastic Column. 
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SOME CONTRIBUTIONS TO THE WEDGE-WATER ENTRY PROBLEM 
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(Proc. Paper 1214) 


SYNOPSIS 


Some of the results of a theoretical and experimental study of the wedge- 
water entry problem are presented. (1) The study is based upon a new formu- 
lation of the basic equations of hydrodynamics especially applicable to this 
problem. The development is given first for the compressible case. This is 
then specialized to the incompressible (water) case and a study of the result- 
ing equations leads to various relations and properties governing the general 
straight-sided wedge water entry phenomenon. 
The new theory is then applied to give the following: 


(a) An exact relaxation solution for an unsymmetrical entry case. 

(b) A closed form analytical expression for the approximate shape of the 
free surface for a thin wedge. 

(c) An approximate solution for blunt-nosed wedges. 


In addition, the theoretical results of (a) and (b) are compared with experi- 
mental results. 


NOTATION 


dimension parameters for blunt-nosed bodies 


direction of normal 

pressure 

(u, v), physical (actual) velocity 

(3,9) = (F ; z ), distance in the conical (3% ) plane 
subscript, “steady state” 


2 210 


Note: Discussion open until September 1, 1957. Paper 1214 is part of the copyrighted 
Journal of the Engineering Mechanics Division of the American Society of Civil 
Engineers, Vol. 83, No. EM 2, April, 1957. 
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time 


3 u,v components of the physical (actual) velocity 

F x,y physical (actual) coordinates of a point 

% F.S. free streamline, free surface 

3 Q q - r, “reduced” velocity 

z R = , a term introduced in a transformation of the flow field 


distance along free streamline, entropy 
a U constant blunt-nosed body velocity 
q Vo constant water velocity for wedge immersion (wedge stationary) 
- a line, a curve in the (3, A ) plane parallel to Q, hence a particle path 


B line, a curve in the ( q ) plane along which £ +s = constant 
3 y 7 conical coordinates 

half wedge angle 

o + id, a plane of transformation 


6 

u 

potential function 
Y stream function 
p 

8 

wW 


density of water 


Q? -i6 , a term introduced for thin wedge, free streamline analysis 
angle of Q line 


gu _ gv 
aw 


( 3 ) Sy the vector operator “del” 


2 


), the physical vorticity 


) the Laplacian operator 


the “reduced” vorticity 


INTRODUCTION 


The wedge-water entry problem was studied first by Wagner(2) who es- 
tablished the similarity property and obtained an approximate solution based 
upon an “expanding plate” analysis. He also gave one value for a total force 
which, presumably, was obtained by an iterative method but no details were 
given. Since then, with one exception, all work on the wedge-entry problem 
has essentially been either an extension or filling-in of Wagner’s original 
presentation (Pierson, (3) Cooper,(4) Karzas,(5) Bisplinghoff and Doherty, (6) 
Schiffman and Spencer(7)), The exception is Garabedian, (8) who utilizes a 
set of Kirchhoff type transformations to obtain a solution for an oblique 
entry subject to unusual initial and boundary conditions. 

In this report, a procedure is used which differs from all of the above. It 
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leads to a new set of equations governing the phenomena. A study of these 
gives various relations and properties for the wedge-entry problem and 

these are presented in this report. Following this, various applications of 
the new theory are given and an experimental verification of some of these 
results is also described. 


The Two Dimensional Unsteady, Non-Viscous, Non-Heat 
Conducting Flow in Conical Coordinates 


It will be shown that the non-viscous, non-heat conducting unsteady flow 
equations, in two dimensions, can be given in terms of new variables, 


7 % 

The flow equations in this form do not contain the length dimensions alone 
and in the ( 3 ) field the flow is represented by a single map which 
represents the flow at all times. This implies that such time dependent flow 
fields grow uniformly with time. If, in addition, the boundary does not 
specifically contain a length parameter, then the field may be called “conical” 
and the entire phenomenon may be called one of “unsteady conical flow.” 
Examples of flows which satisfy all of the above conditions are: 


(a) A shock striking an infinite straight sided wedge. 

(b) An infinite straight sided wedge travelling with constant velocity and 
entering an infinite ocean, or conversely, an infinite ocean travelling 
with constant velocity and engulfing an infinite straight-sided wedge. 


The case (a) will be discussed first. This corresponds to compressible 
flow conditions. The equations will then be specialized to incompressible 
flows and will therefore be the equations governing the phenomena of type (b). 

Consider the two dimensional, unsteady, non-heat conducting non-viscous 
compressible flow, with 


If it is assumed that 2 and U~ can be given in terms of the new variables 


(3 ? ) defined by 
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and it will be found that the equations are independent of x, y and t. Because 
the boundary is also independent of any 7S parameter, it follows that the 
flow pattern can be represented in the (2 » plane. 

The basic equations are 


conservation of mass 
A "ty +PUx +P Vy = 
conservation of momentum 
equation of state 


(PS) (3) 


and for changes of state which are adiabatic, an equation of conservation of 
energy, in terms of entropy S, 


7% Sy =O (4) 


Hence, corresponding to five unknowns, u, v, p, p, and S there are five 
independent equations. 
Now introducing the coordinates ( 3 ” ) and noting that 


and defining a “reduced velocity,” Q, by means of 


Q= 
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the five equations become 


conservation of mass, )+ =o (7) 


conservation of momentum, Q:7Q + Q = 3 Ve (8) 2 


state, = CF, S) (9) 


conservation of energy, Q- VS=o (10) 
OF» 
‘ Equations (7) through (10) are the general compressible flow form of the . 
= basic equations in terms of the reduced velocity Q, and the variables 7 


3, y ). 


Py Some Properties of the Wedge Water Entry Flow Field 
we in Conical Coordinates 


For the incompressible flow, with which this report shall hereafter be ost 
om concerned, the above five equations specialize to three equations in terms of 4 
“am 6h uu, v, and p. The first of these, the conservation of mass, leads to 

6. I The conservation of mass equations becomes, in terms of Q and ( 3,*) ), 


7-Q (11) 


This follows directly from equation (7) with constant p . 


; The second and third of the basic incompressible equations are the conser- 
J vation of momentum equations. Eq. (8) becomes, for constant density 
CPQ-VQ +PQ --vPp (12) 
Using the familiar vector identity 
z 
)- @x(7xQ) (13) 
Eis equation (12) becomes, upon collecting terms 
For an irrotational motion WX Q zO . Hence 
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q II—The conservation of momentum equation for the incompressible fluid in 
= irrotational motion is given by 


(15) 


A stream function, Ww , may be defined as follows 
(16) 


which identically satisfies Eq. (11) as may easily be verified. Taking VX 
of Eq. (14) and using Eq. (16), there follows 


I1I—The incompressible fluid in rotational motion satisfies the equation 


ig 
a This is equivalent to 
VV ¥)-7¥=0 
Also, if we define a “reduced vorticity” by (Z= 7X Q then 
=z 
V¥ =22 
4 and it follows from Eq. (18) that 
. ‘% IV—The incompressible fluid in rotational motion also satisfies the equation 
d Q- =O (20) 
4 Note that 2 is equal to the physical vorticity, Ww , multiplied by t. 
Also, 
a V-—If an incompressible, conical fluid motion has Q at any time, the fluid 
z motion is always irrotational. This is the Helmholz vorticity statement 
‘ for the incompressible, conical flow. 
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If the fluid motion is irrotational, then it follows from Eq. (16) that 


ViI—The stream function, , Satisfies the Laplace equation in ( 3 ’ 7 ) 
coordinates everywhere in the field, or 


7’°¥=O (21) 


and Eq. (17) is identically satisfied. 

In the present study, the fluid motion is irrotational before wedge entry. 
Therefore it is also irrotational after wedge entry. Because the motion is 
irrotational, there exists a potential, y , Such that 


Eq. (15) can also be stated in the form 
Q* 
VCS+E 
Therefore 


and from Eq. (22) and because of irrotationality 


7°9=0 (26) 


Also, in the ( 3 by ) plane, lines of constant f are perpendicular to lines 
of constant ¥ . All of the above is stated as 


VII—In the (3 v ) plane there exists a potential function defined by 


oF — 
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such that 


(26) 


Also, lines of f = constant are normal to lines of ¥ = constant. 

In analogy to streamlines of steady flow we can introduce, in the ( 3, ”») 
plane, the concept of Q lines. These lines, which are everywhere tangent to 
the direction of Q, may be defined as a family of lines 


Xx(Z ” ) = constant (27) 
This being so, it follows from Eq. (15) that 
VIlI—In the ( 3, 7 ) plane, if we define 

- 
then there are lines 

AC ») = constant (29) 
which are everywhere perpendicular to the lines 

ox = constant (27) 

Thus, up to this point, two mutually perpendicular sets of curves in the 

( Ey ) plane have been defined for the general non-viscous, non-heat con- 
ducting, incompressible fluid in irrotational motion. These are 


POTENTIAL CURVES § = CONSTANT 


“STREAMLINES” 7 = CONSTANT 


= CONSTANT 


6 (4,9) = CONSTANT 
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Also, up to this point the equations and relations obtained are perfectly 
general and held without regard to origin of coordinates. The discussion will 
now be tied into the particular problem being considered. The major portion 
of this report is concerned with the behavior of an infinite straight sided 
wedge travelling with constant vertical velocity Vo, and entering an infinite 
ocean with undisturbed surface horizontal and extending to infinity in both di- 
rections. The pressure in the air and in the ocean before immersion is 
zero. 

For our purposes, it will be convenient to assume the problem in the fol- 
lowing equivalent form (see Fig. 1). The wedge is stationary, the infinite 
ocean is moving vertically upward with a constant velocity V, and engulfs 
the wedge. The origin of the ( 3, ) coordinates will be taken at the point 
of the wedge. After engulfing the wedge, the ocean has free surfaces AB 
and A'B' which attach to the wedge at points A,A'’ as shown, which are at 
zero pressure and which approach the undisturbed surface position at infinity 
tangentially. The velocity and pressure of the ocean at infinity are Vo and 
zero. Along the wedge face, the fluid velocity is tangent to the wedge surface. 
Hence the wedge face is a Q line. The motion is irrotational. It may also be 
noted that continuity of mass requires that the volume of fluid displaced 
above the initial position of the free surface must equal the immersed volume 
of the wedge. Unless otherwise stated all of the following developments refer 
to the conical ( 3, ” ) plane. 


IX—From Eq. (15) it follows that a Bernoulli type equation in the ( 3 
plane for the flows being considered, is (S being the distance LE RP 
along Q lines), 


(30) 


This is an expression for the variation of Q’*/2 + p/p alonga Q line. 
The corresponding steady state Bernoulli equation is (S being the distance 
measured along streamlines in the physical (x,y) plane), 


+#) 


The following is an important relation: 


X—In the ( Z ) plane, Q lines are particle paths. 


__ The proof follows: The particle path as defined in the physical plane, 
R(x,y), with q(u,v) is given by 
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Thus, in the ( J, ” ) plane, Q lines are the particle paths. It is interesting 
to note that if, in addition, 


CT: log t 


Q 


so that in the time distorted = plane, the Q lines are not only the particle 
paths, but the particles move along them with velocity Q. 

The three statements which follow immediately have been proven elsewhere 
(as noted) and are included herein because they are utilized in proving some 
of the relations which come after. 

Several proofs have been given for XI. See, for example, Garabedian. (8) 
XII was proved by Cooper(4) and XIII is a classical result given, for example, 
in Lamb. 


XI—The distance between particles on the free surface, in the (x,y) plane re- 
mains constant with time. 


XII—The free surface particle in contact with the wedge face at the initial in- 
stant of immersion is always in contact with the wedge face. 


XIlI—For the flow conditions considered herein, a particle which is in the 
free surface at any time is always in the free surface. Hence, the free 
surface is a particle path, or a Q line, and we shall speak of it as the 
“free streamline,” i.e., F.S. 


It can now be shown that 


XIV—The free streamline in the ( Zz » ) plane attaches to the wedge at a 
finite distance from the wedge point. 


The proof is as follows: 


(1) The wedge is stationary and the infinite ocean is travelling upward with 
uniform velocity Vo at infinity. 

(2) The conical property requires that the complete flow field be established 
at the instant after immersion, t = 0°, and that this field be similar to 
itself at all times thereafter. 
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(3) According to XI, the distance between any two particles on the free 
streamline is constant, starting from the time t = 0 when the wedge 
begins immersing. Also, according to XII, a particle in the free 
streamline at any time is always in the free streamline. 

(4) If the free streamline does not attach to the wedge at a finite distance 
from the wedge point (i.e., attaches at infinity), then because of (a) 
above, at least one particle on the free streamline must have had an 
infinite velocity at the instant of immersion t = 0. 

(5) If one particle had an infinite velocity, then because of (3) all particles 
in the free streamline have an infinite velocity at the instant of immer- 
sion and therefore at all times. This violates the boundary condition 
in (1) and is therefore not possible. Hence the free streamline must 
attach to the wedge at a finite distance from the wedge point. 


From XIV and because of XII and the property of the conical flow plane, 
it follows that 


XV—The free streamline particle at the wedge attachment point, A, travels 
with a constant velocity q,, which is equal to r,, the distance of this 


particle from the wedge point in the ( 4,7) ). Therefore, for this 
particle Q, = 0. 


It is now possible to derive the boundary condition which must be satisfied 
by Q along the free streamline. The complete statement is 


XVI—At any point on the free streamline in the ( Z y) ) plane 


(32) 


s,| = Sp 5, 


where |Q, s | is the value of Q at this point on the free streamline, and 
Sp s is the distance, measured along the free streamline, from the wedge 
attachment point to the point in question. In addition, of course, the vector Q 
is everywhere tangent to the free streamline. 

The last part of XVI follows at once because of XIII. The proof for the 
first part follows—see Fig. 2. 


(1) At point A, 
(2) At point B, r = rp as shown 


(3) Between A and B, 4q must be perpendicular to Ar as shown in the 
figure. This is so because the pressure gradient on the free stream- 
line is perpendicular to the free streamline 


(4) Hence qp 4 Tp since from (1) above, qp will equal Tp only if 
= ( Map 
(5) Hence Qp #0 
(6) On the free streamline, p = 0, and therefore from Eq. (15) 
Q(dQ/dS + 1) = 0 


everywhere along the free streamline, where S is in the direction 
tangent to the free streamline. 
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(7) At B, Qp # 0, so that from (6) above, | Qp| 2 | 4 r | and continuing 
along the free streamline to any point on the free streamline, we obtain 


= Sps 
A result of some interest is 


XV0—Throughout the ( 3 9) plane flow field the function (3 1 ) 
satisfies the equation 


Vp 

This follows from 
I 
and 7/3 


Take V .- (15) and using Eq. (11), Eq. (33) is obtained at once. 
Also, we may show that 


XVI0-—In the (3, ” ) plane, along the wedge face, the pressure will be a 
maximum or a minimum only if 


Q:0 
Proof: From II 
Q= 
$+ ) +Q=0 
and because the wedge face is a Q line, this becomes, along the wedge 


+Qd5=0 


=P 4+QdQ #Qd5=0 


SP +Q (dards) =o 


From this 
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(37) 


But, along the wedge 


An Exact Relaxation Solution for an Unsymmetrical Wedge 


In view of equations (22), (25) and (26), it is possible to set up a numerical 
relaxation procedure for obtaining exact solutions to wedge-water entry 


problems. 
The solution requires finding a which satisfies 


ae a = 0, everywhere in the field and 
the following boundary conditions: 
f - Q’/2 - r*/2 along the free streamline 


(It can be noted at this point that the free streamline at infinity has a value 


i 
for ff given by a 5 ) 


of = 0 along the wedge face, n being the direction 
on ~ of the normal to the wedge 


af = - v, at infinity 


The problem is complicated by the fact that the free streamline is not 
known beforehand. In this respect the difficulty is similar to the one which 
appears in the problem of percolation of a fluid through a porous medium, 
which was solved by relaxation methods by Shaw and Southwell. 10) In es- 
sence, the procedure in these cases is to assume a free surface and to obtain 
a solution which satisfies the field equation and the boundary conditions sub- 
ject to the assumed free surface. Then, using this known solution, there will 
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be another equation or relation which leads to an independent check of the as- 
sumed surface shape. In general, the assumed shape will not check out, so 
that another shape must be assumed and the process repeated. 

In the present problem, the procedure formalizes to the following: 


(1) Assume a free streamline shape. 
(2) Calculate f along the free streamline, using 


= Q’/2 - r’/2 


(3) Obtain a relaxation solution of the Laplace equation subject to the given 
boundary conditions. 

(4) Having obtained everywhere, it is now possible to determine the 
physical velocity,’ q, along the free streamline. Hence, at all points on 
the free streamline, Q = q - TF can be obtained. These values of Q 
must check the actual assumed values of Q. If they do not check, then 
the free streamline shape must be altered and the entire process re- 
peated. 

(5) Having the correct solution for f, it is now possible to obtain values 
of p/p along the wedge face, since 


= Q’°/2 + p/p - 
and Q and r are known. 


This method is now applied to the unsymmetrical entry of a 90° wedge as 
shown in Fig. 3. Because the Laplacian is invariant with respect to rotation 
of axes and because the relaxation numerical work is simplified thereby, the 
net lines are taken parallel and perpendicular to the wedge face as shown. 

As initial guides for the assumed free streamline, the two continuity con- 
ditions on displaced volume and free streamline length are used. In the be- 
ginning, rather coarse nets are used. These permit checks on the F.S. and 
they are used until an acceptable check is obtained. Then the net is made 
finer and the process repeated. 

In one of the early solutions, it was found that the fluid along the wedge 
face OA, in the neighborhood of 0 is moving in the direction of 0. This im- 
plied the existence of an infinite velocity singularity at 0 corresponding to 
flow about a sharp corner. Methods of handling this are described by 
Motz‘11) and these were used. In the present problem, it was found neces- 
sary to use a fine net in the neighborhood of the sharp bend of the free 
streamline. The final solution is shown in Fig. 3, which shows 


(a) The distribution of P/5 p Vs along each wedge face. 


(b) The free streamlines 
(c) The network used in the relaxation solution. 


An experiment was performed to check the results. A wedge with plastic 
end plates (to simulate two-dimensional effects) was dropped into a tank of 
water at constant velocity and a high speed flash photograph was taken of the 
resulting disturbance. For the case under consideration, the photograph is 
shown in Fig. 4. 

In connection with this photograph, it must be remembered that the theory 
assumed a perfect fluid in which gravity effects were neglected. For the 
most part these assumptions are reasonable for the experiment performed. 
However, in the neighborhood of the free streamline-wedge attachment point, 
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it is obvious that gravity has caused “fall back” of the particles. It is, how- 
ever, still possible to obtain a comparison between the theoretical shape of 
the free streamline and the experimental shape, if the photographic free 
streamline is simply corrected for its gravity trajectory fall back at its 
outer end to give the wedge attachment point. This has been done, and the 
resulting free streamlines are shown in Fig. 5. 

It will be recalled that the theoretical solution called for an infinite veloci- 
ty (and hence negative infinite pressure) at the wedge point. These are, of 
course, impossible to obtain experimentally. Thus, some disagreement be- 
tween theory and experiment is to be expected on this score. 

Also, the photograph Fig. 4 clearly shows leakage of fluid between the 
plastic end plate and the wedge on the high pressure side of the wedge. This 
would reduce the amount of fluid on this side of the wedge and the comparison 
between theoretical and experimental free streamlines at this face clearly 
shows this effect. 

On the whole, however, the agreement between theory and experiment for 
this unsymmetrical wedge is good. 


The Free Streamline Shape for the Thin Wedge 


The theory of this study may be utilized in obtaining an approximate solu- 
tion for the thin wedge free streamline shape. The method is analogous to 
Kirchhoff’s free streamline theory of steady incompressible flows. The so- 
lution which will be obtained satisfies the following conditions: 


(1) The boundary condition on Q along the free streamline 

(2) The F.S. will be tangent to the wedge at its attachment point. 

(3) The free streamline asymptotic condition, i.e., the physical require- 
ment that at infinity the free streamline approaches the undisturbed 
surface level asymptotically. 

(4) The volume displacement condition, i.e., the requirement that the im- 
mersed volume of the wedge must be equal to the volume of the fluid 
displaced above the initial horizontal water surface. 

(5) Because in the limit as the wedge angle becomes zero the radius of 
curvature at the free streamline wedge attachment point is zero, the 
solution which will be obtained is one for which the radius of curvature 
will be zero at this point. 


Let (see Fig. 6) 
6 the half angle of a symmetrical thin wedge 
~ inclination of the Q lines with respect to vertical axis. 


Introduce the quantity 


Consider the term 
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in which S is the distance measured along each Q line from the point where 
each line intersects the wedge. 
Then 
log R = log Q + i0 (42) 


Along the free streamline, S = Q exactly, and therefore along the free 
streamline 


log R = i6 (43) 


The actual flow field of Fig. 6 will now be approximated on a log R plane 
as shown in Fig. 7. 

AB represents the free streamline exactly as f— 0. BEo represents 
the flow at infinity between B and some point Eo such that between these 
points 

(1) 6 = 7/2 and 

(2) S/Q goes from unity to oo between B and Eo. Between E and 

D oo it is assumed that 
(1) @ changes from 7/2 to 0 and 
(2) S/Q remains infinite. 


The entire half flow field thus lies within the rectangle Eno BADoo in the 
approximation assumed here. 


Transform the log R plane of Fig. 7 intoa w = o + id plane, Fig. 8, 
such that E,, BAD, lies on the real axis, A = 0, and so that B is at 
o@ = 0 and A is at @ = +1. The transformation that accomplishes this is 
log R = K cosh™* w + L (44) 
in which K and L are constants which determine the origin and the scale of 
the transformation. Using the mw plane values for the locations of A and B, 
this becomes 
log R = cosh™* (45) 
Also, along AB, log R = i@, so that along AB this becomes 
= cosh’ (46) 
= cos (47) 


du=-sin@da6 (48) 


(49) 


Also, along AB d3Z = cos @ ds 
sin 6 dS 


In order to obtain an equation for the free streamline AB which is con- 
sistent with all the above, it is necessary that some relation be established 
between S and /. In particular, the relation will be given in terms of 
dS/du. 
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First note, 


and from (47) 
d6@/du = - 1/sin 


Also, it is necessary that 


0 at A, in order that the radius 
of the curvature be zero here 
dS/d@ = 
co at B, in order to satisfy the 
asymptotic condition 
All these are satisfied by taking for dS/d# the simple relation 
dS/du# = M/p = M/cos 8 (53) 


where M is some constant which will be determined in order that the con- 

stancy of volume condition be satisfied. See Fig. 9. Note also from Eq. (57) 

which follows that M represents the maximum rise of the free surface. 
From (49) through (53) 


M sin @ d 8 
Also, the displaced half-wedge volume is V6 6/2. 


Now equating the volume under the free streamline and above the undis- 
turbed free streamline to the displaced half-wedge volume, one obtains 


fom v2 $2 


or, substituting the values from (54) and integrating, 
M = - Vo V2 (56) 


Integrating the expressions for d 3 and a the parametric equations of 
the free streamline are finally obtained as 
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- Vo V20/n (1 - cos 8) 
(57) 
7m (1/2 10g ++ 8 gin @) 


1 - sin 8 


in which and are here measured from the origin as shown in Fig. 9. 

An experiment was performed to obtain a check on this result. A thin 
wooden wedge, of total angle 20°, with plastic side plates, was dropped into a 
tank of water with constant velocity Vo. Photographs were taken of the re- 
sulting disturbance. In Fig. 10, the experimental free streamline for this 
wedge angle is compared to the theoretical free streamline given by Eq. (57). 

In connection with this figure, it must again be recalled that the basic 
theory of this investigation assumed a perfect fluid—i.e., a non-viscous, 
gravity free, surface tension free fluid. For the most part, these assumptions 
are reasonable. However, in the neighborhood of the free streamline attach- 
ment point it would be expected that viscosity and surface tension effects 
might be important. In the present case, the experimental free streamline 
had a reversal of curvature near the attachment point—which is impossible 
in a perfect fluid and is probably due to surface tension effects. 

On the whole, however, the agreement between the theoretical and experi- 
mental free streamline shapes is good. 


Blunt-Nosed Wedge Water Entry 


Because of equations (22), (25) and (26) it is possible to obtain approximate 
solutions for straight-sided wedge entry problems as follows: 
We require a 3,” ) function which satisfies 


(1) v4 = 0, everywhere in the field 


everywhere in the field 


2 2 
(3) § = om - 7 along the F.S. 


2 
(4) g oi Yo at infinity on the F.S. 
(5) dn = O along the wedge face. 
(@) 23 =- Vo at infinity. 


§ functions which satisfy all of the above except (3) may be obtained by 
utilizing the known results of steady-state potential flow theory of classical 
hydrodynamics. Having everywhere in the flow field, we can use equation 
(25) to find the pressure along the wedge face or any other desired property. 

This method will give rapid, approximate solutions for straight-sided 
wedges such as were considered in the previous portions of this paper. We 
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need only use the known potential solutions for steady flow about diamond 
shape wedges. 

It is possible to apply this method to obtain approximate solutions for un- 
steady water entry flows for two dimensional blunt-nosed bodies, for which 
the classical theory of potential flow will permit the development of prac- 
tically any desired shape. (12) Typical results obtained using this approxima- 
tion are: 


(1) Fig. 11 shows the pressure distribution along a blunt-nosed wedge 
shape developed by a vortex pair in a uniform flow field.(13) The 
region of negative pressure will exist only in an ideal fluid. We may 
consider it as a cavitation or separated region for the real fluid. 

(2) A blunt-nosed body may be developed by a steady potential flow source- 
sink pair in a uniform field. 14) One half of this immersing symmetri- 
cal body could be the nose of a two dimensional projectile or similar 
body. When this nose is just fully immersed, there is an unsteady state 
pressure, p, at the nose tip which is greater than the steady potential 
flow stagnation pressure, ps. The amount of this excess is given by 


b (b + a) 
i ae log ( ) (58) 
5p U ab (b - a) 


If a/b — 0, then the nose shape becomes a semi-circle and the expres- 
sion becomes 


P- Ps 
pu’ 


i = 2 (59) 
2 


The same method can be used to give approximate solutions for unsym- 
metrical two dimensional blunt-nosed body entry and for three dimensional 
shapes as well. An extension of the method permits one to obtain approxi- 
mate solutions when the entry velocity is variable. 


CONCLUSION 


The paper presented some of the results of a study of the wedge-water 
entry problem based upon a new formulation of the basic equations of hydro- 
dynamics applicable to this problem. Included were 


(a) a development of the equations and the various relations and properties 
which follow from these. 

(b) an exact relaxation solution for an unsymmetrical wedge. 

(c) a closed form analytical solution for the approximate shape of the free 
streamline for a thin wedge. 

(d) an approximate solution for blunt-nosed wedges. 

(e) comparisons between theory and experiment for cases (b) and (c). 
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Fig. 4 High speed photecraph ef the unsymmetrical 
entry ef the 90° wedge 
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Fig. 6 The thin wedge 
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Fig. 7 The log R plane in thin wedge analysis 


4 

4 % 

7) 

A a 

ian 

FS, 

Log R Fame. 

4 

A <2 = 

me. 

A Deo a 

1 @ Fig 8 The r plane in thin wedge analysis Ae 

4 


: ASCE BORG 1214-27 


{ Fig. 9 The free streamline in thin wedge analysis 
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Fig. 10 Comparison of theoretical and experimental 


free streamlines for the thin (20°) wedge i 
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Discussion of 
“THE VISCOUS SUBLAYER ALONG A SMOOTH BOUNDARY” 


by H. A. Einstein and Huon Li 
(Proc. Paper 945) 


H. A. EINSTEIN, ! M. ASCE, and HUON LI.—Both discussers, Messrs. E. 
Silberman and N. N. Ambraseys have helped the writers in one vital point. 
They made it very clear that the problem is difficult to handle experimentally 
and theoretically alike. The experimental difficulty is caused predominantly 
by the turbulence outside the sublayer which moves past any place of obser- 
vation at high velocity and causes pressure fluctuations of equal magnitude 
but with higher frequency then these developed in the sublayer itself. The 
theoretical difficulty is caused by the fact that the sublayer should constitute 
a transition between the still boundary and the turbulent velocity profiles 
while a mathematical solution for the transition could be found only between 
the still boundary and a uniform flow of constant velocity. 

No objections were voiced against the main ideas of the paper: 

1) that the sublayer was interpreted as a non-steady flow with an average 
flow pattern rather than as a steady flow with superimposed fluctuations. 

2) that in areas of predominantly viscous character turbulent effects may 
be neglected (build-up phases) and that in areas of turbulent character the 
viscous effects are negligible (mixing with main flow). 

These two ideas appear to be basic and new and may prove to be helpful in 
the description of other similar phenomena, too. 

Mr. Silberman proposes the use of another definition for the effective sub- 
layer thickness. He is welcome to change that definition at will since the 
choice is completely arbitrary. The proposed use of dye was actually applied 
by the authors to show that at some instances a distinctly laminar flow pat- 
tern exists near the wall while at some other instances turbulence appears to 
reach all the way to the wall. The processes are too fast for visual observa- 
tion but show up very well in high speed photography. 

Mr. Ambraseys’ observations of the statistical flow description in the 
range of the sublayer seem to be most interesting. Unfortunately, no obser- 
vations are available from the experiments described in this paper that would 
permit to check Mr. Ambraseys’ results. 

In conclusion, the writers like to extend the heartiest thanks to the dis- 
cussers for their help and suggestions and they hope that more measurements 
may become available soon that will permit a more detailed description of 
the sublayer with the help of the proposed theory. 


1. Associate Prof. of Mechanical Eng., Univ. of California, Berkeley, Calif. 
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Discussion of 
“LATERAL BRACING FORCES ON BEAMS AND COLUMNS” 


by William Zuk 
(Proc. Paper 1032) 


CORRECTIONS.—Figs. A and B to which Prof. Posey referred in his dis- 
cussion of this paper were omitted from Proc. Paper 1154. These figures 
are reproduced herewith: 
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WILLIAM ZUK, | .—The writer wishes to thank Messrs. Hauer, Posey, and 
Segedin for their interest in commenting on the subject of lateral bracing. 

Both Messrs. Hauer and Posey noted a lack of discussion on the stiffness 
factor in connection with elastic bracing. The writer is quite aware of the 
role of stiffness and has studied the subject at length, even though it was only 
briefly mentioned in the paper. Stiffness was deliberately de-emphasized in 
the paper because, as Mr. Hauer himself pointed out, “the stiffness require- 
ment for the bracing member is not hard to fulfill.” This conclusion was also 
reached by the writer for many other practical examples. It was felt there- 
fore that a force analysis would more adequately serve the needs of the prac- 
ticing engineer. Stiffness, of course, has no meaning for cases of immovable 
supports, where forces themselves must be had. 

Another comment of Mr. Hauer concerning the doubtfulness of using a 
crookedness parameter calls for reply. To those familiar with the buckling 
phenomena, it is clear that the initial imperfection is one of the most funda- 
mental parameter. As uncertain as this crookedness may be, it neverthe- 
less—in an analytical study—must be presented as the parameter it is. 
Reference is made to the classic “Secant Column Formula” for column 
buckling which also has in it an imperfection parameter. The “Secant 
Formula” stands as one of our most accurate column formulas. 

Unfortunately, Mr. Posey’s drawings were accidentally omitted from his 
discussion but are shown just prior to this closure. Nevertheless, the writer 
would like to make a side note that in a continuous column with a brace at 
mid height, the true position of the flex point is not exactly at the center. 

Mr. Posey’s analysis however gives a very simple and effective picture of 
column-bracing action. 

To Mr. Segedin, the writer is endebted for his refinement of the solution 
of Case I and to his very interesting extension to the case of moment bracing 
of a column. 

In closing, the writer wishes to emphasize the importance of the portion of 
the paper on beam bracing which is the more difficult subject and on which 
experimental work should be done. 


1, Associate Prof. of Civ. Eng., Univ. of Virginia, Charlottesville, Va. 
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